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ABSTRACT 


The  objective  of  this  project  was  to  carry  out  both  experimental  and  analyt¬ 
ical  studies  leading  to  improved  design  procedures  for  predicting  thrust  termination 
of  single -chamber  controllable  solid  motors.  The  experimental  work  consisted  of 
measuring  quantitatively  the  effects  of  incident  thermal  radiation  on  low  pressure 
burning  rates  and  deflagration  limits  of  typical  solid  propellants .  A  technique  was 
developed  in  which  small  cylindrical  samples  of  the  propellant  were  burned  inside 
an  electrically  heated  tube  furnace.  At  furnace  wall  temperatures  up  to  1500°F, 
the  burning  rates  were  observed  to  be  as  much  as  50%  greater  at  lOpsia  than  they 
were  when  surrounded  by  room  temperature  walls.  Analytical  work  was  performed 
using  an  improved  mathematical  model  of  the  transient  combustion  process  of  a 
solid  propellant.  This  study  led  to  the  conclusion  that  experimentally  characteriz¬ 
ing  the  extinguishability  of  a  propellant  in  terms  of  a  critical  dp/dt  provided  very 
little  useful  design  information.  On  the  other  hand,  characterizing  extinguishability 
in  terms  of  the  product  L^r^  was  shown  to  provide  a  design  criteria  with  the  prom¬ 
ise  of  very  general  applicability. 
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INTRODUCTION 

The  accurate  prediction  of  the  thrust  of  a  rocket  motor  during  the  transient 
periods  that  precede  or  follow  steady-state  operation  is  a  difficult  problem.  Even 
with  simplifying  assumptions,  the  mathematics  describing  the  conservation  of  mass 
and  energy  are  complex  and,  more  importantly,  the  burning  characteristics  of  the 
propellant  under  transient  conditions  are  not  fully  understood.  Difficulty  notwith¬ 
standing,  obtaining  satisfactory  solutions  to  this  problem  is  becoming  more  critical. 
This  aspect  of  design  is  especially  important  for  controllable  motors  or  for  conven¬ 
tional  motors  having  thrust  termination  capabilities .  The  need  for  improved  design 
capabilities  is  indicated  most  forcefully  by  the  recent  failures  of  scaled-up  research 
motors  to  completely  terminate  thrust  on  command  (1).  Instead  of  shutting  down  as 
had  been  predicted  on  the  basis  of  small-scale  data,  the  motors  continued  to  burn  at 
a  low  level. 

The  research  reported  herein  was  designed  to  provide  information  to  improve 
the  ability  of  the  motor  designer  to  accurately  predict  marginal  extinguishment  con¬ 
ditions.  The  work  was  planned  around  two  important  observations:  (1)  the  steady- 
state  low  pressure  strand  burning-rate  data  that  is  often  used  for  motor  design 
purposes  can  be  seriously  in  error  because  of  thermal  radiation  effects,  and  (2) 
correlations  can  be  made  which  indicate  that  the  extinguishability  of  a  wide  variety 
of  propellants  increases  inversely  with  the  steady-state  burning  rate.  The  first 
observation  indicated  a  need  for  improved  methods  for  obtaining  laboratory  data  . 

The  second  observation  indicated  a  need  for  analytical  work  leading  to  more  accept¬ 
able  general  correlations. 

The  work  accomplished  can  be  logically  divided  into  four  phases:  Experi¬ 
mental  Work,  Data  Compilation  and  Correlation,  Parametric  Design  Studies,  and 
Recommended  Design  Procedures. 

It  should  be  noted  that  the  analytical  work  was  directed  specifically  toward 
applications  in  the  design  of  single -chamber  controllable  motors;  however,  the  con¬ 
clusions  regarding  low  pressure  burning  rates  and  the  sizing  of  the  nozzle  area  can 
also  be  applied  in  the  design  of  the  aft  chamber  of  a  dual-chamber  controllable  motor. 
The  problem  of  spontaneous  reignition  following  thrust  termination  was  not  consid¬ 
ered  during  this  study . 
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II 


EXPERIMENTAL  WORK 


1 .  Objectives 

One  design  technique  employed  to  insure  extinguishment  of  single -chamber 
controllable  solid-motors  is  to  size  the  nozzle  throat  area  such  that  with  the  nozzle 
fully  open  the  motor  L*  and  pressure  will  both  be  less  than  some  critical  value. 

The  magnitudes  of  these  critical  values  are  determined  empirically  using  test  fir¬ 
ings  of  small  motors  (2). 

Problems  have  been  encountered  in  the  use  of  this  technique  because  the 
steady-state  burning  rates  determined  in  small  motors  or  with  strand  burners, 
which  are  used  to  compute  the  full-scale  nozzle  throat  area,  can  be  significantly 
less  than  the  burning  rates  that  occur  in  the  full-size  motor  (1).  This  discrepancy 
in  burning  rate  can  lead  to  continued  burning  after  the  nozzle  pintle  has  been  with¬ 
drawn,  rather  than  extinguishment. 

The  objective  of  this  phase  of  the  program  was  to  investigate  methods  for 
reliably  characterizing,  in  the  laboratory,  the  low-pressure  burning  rates  and 
deflagration  limits  that  might  be  expected  in  full-scale  motors.  It  was  reasoned 
that  the  principle  difference  between  the  combustion  process  in  full-scale  motors  and 
in  small  motors  or  strand  burners  was  the  net  supply  of  thermal  radiation  to  the 
burning  surface  of  die  solid  propellant.  In  the  full-scale  motor,  the  surface  would 
be  exposed  either  to  additional  burning  surface  or  to  the  hot  insulation  covering  the 
nozzle  and  pintle  housing.  Thus  the  net  supply  of  thermal  radiation  incident  to  the 
burning  surface  would  be  near  zero,  or  positive  if  the  insulation  surface  were  at  a 
higher  temperature  than  the  propellant  surface.  On  the  other  hand,  in  the  strand 
burner  the  burning  surface  would  be  exposed  to  the  cold  walls  of  the  pressurizing 
container  and  the  net  supply  of  thermal  radiation  would  be  negative.  In  small  motors 
the  ratio  of  burning  surface  area  to  the  total  surface  area  of  the  combustion  chamber 
is  usually  considerably  less  than  it  is  in  full-scale  motors,  and  the  exposed  area  is 
not  so  well  insulated.  Thus,  the  net  incident  thermal  radiation  to  the  burning  sur¬ 
face  may  also  be  negative  in  small  motors. 


2.  Apparatus 

The  temperature  at  the  surface  of  burning  propellants  has  been  measured 
by  various  means  (3,  4)  and  typically  has  been  reported  to  be  near  1000-1100°F 
(550-600°C).  Therefore,  if  the  large  motor  environment  is  to  be  simulated  in  a 
laboratory  strand  burner,  the  strand  should  be  exposed  to  inert  surfaces  heated  to 
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temperatures  in  this  range,  or  higher,  or  the  propellant  sample  should  be  prepared 
in  such  a  way  that  the  burning  surface  is  exposed  principly  to  additional  burning  sur¬ 
face  . 


In  exploratory  experiments  conducted  at  BYU  prior  to  this  program,  strands 
of  an  ammonium  perchlorate -oxidized  composite  propellant  were  burned  both  in  the 
conventional  manner  and  also  by  placing  two  strands  end  to  end  so  that  their  burning 
surfaces  were  opposed.  At  1  psia  the  burning  rate  of  single  strands  was  0.028  in/ sec 
and  the  burning  rate  of  the  opposed  strands  was  .053  in/sec .  This  large  difference 
in  burning  rates  could  have  been  partly  due  to  differences  in  convective  heat  transfer 
resulting  from  the  different  velocities  of  the  combustion  products  relative  to  the 
burning  surfaces.  Similar  uncertainties  were  encountered  using  different  sample 
configurations  to  simulate  motor  conditions.  It  was  decided  early  in  the  program, 
therefore,  to  construct  a  strand  burner  with  heated  walls  in  order  that  the  effects 
of  incident  thermal  radiation  might  be  isolated. 

Figure  1  shows  a  diagram  of  the  apparatus  that  was  constructed.  The  basic 
part  of  this  apparatus  is  a  tube-furnace,  Model  No.  423,  manufactured  by  Electro 
Applications,  Incorporated.  This  furnace  is  rated  at  1200  watts  and  with  an  alumina 
tube  can  be  operated  at  temperatures  up  to  2750°F.  In  initial  tests  with  propellant 
samples  inside  the  furnace,  the  alumina  tubes  fractured  due  to  thermal  shock  follow¬ 
ing  ignition  of  the  sample.  Consequently,  the  alumina  was  replaced  with  a  tube 
made  of  stainless  steel.  This  tube,  which  is  2-1/4"  O.D.  with  .065"  wall,  type  347, 
has  performed  satisfactorily.  Because  of  its  greater  thermal  conductivity,  however, 
more  heat  is  conducted  to  the  cooling  coils  at  the  sides  of  the  furnace,  and  the  max¬ 
imum  operating  temperature  is  approximately  1500°F. 


3.  Experimental  Procedure 

Prior  to  placing  the  propellant  sample  into  the  heated  section  of  the  tube, 
this  section  is  allowed  to  reach  the  temperature  specified  for  the  test.  Electrical 
current  is  supplied  to  the  heating  wires  of  the  furnace  through  an  automatic  controller 
which  employs  a  thermocouple  sensor  inside  the  heated  section.  The  controller 
automatically  adjusts  the  electrical  current  to  maintain  the  desired  temperature. 

The  tube -furnace  is  connected  through  a  large  tank  to  a  high  capacity  vacuum 
pump.  The  pressure  in  the  tank  is  regulated  by  controlling  air  leakage  into  the  tank 
through  a  bleed -valve.  A  small  flow  of  nitrogen  is  admitted  to  the  tube  to  provide  a 
continuous  purge  during  the  test. 

The  propellant  sample  was  placed  in  the  heated  section  with  a  hand-operated 
push  rod.  Immediately  after  it  had  reached  the  proper  position,  controlled  by  a  pin 
on  the  push-rod,  it  was  ignited.  Ignition  was  accomplished  with  an  electrically- 
heated  nichrome  wire  in  contact  with  propellant  surface,  the  surface  previously 
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being  coated  with  a  paste  consisting  of  potassium  perchlorate,  ammonium  per¬ 
chlorate,  titanium,  boron,  and  polyisobutylene. 

The  burning  time  of  the  sample  was  detected  with  a  photocell,  Raytheon 
EM1502,  mounted  in  the  end-cap  of  the  tube.  The  photocell  signal,  along  with  the 
signal  from  a  pressure  transducer  also  mounted  in  the  end-cap,  was  recorded 
using  a  Honeywell  Model  1508  Visicorder. 

Several  different  sample  configurations  were  tried.  The  configuration  pro¬ 
viding  the  most  reliable  data  was  a  solid  cylinder,  0.75  inches  in  diameter,  with 
lengths  ranging  from  0.1  to  0.6  inches.  The  ends  of  these  samples  were  machined 
to  insure  that  they  were  parallel  and  that  accurate  measurements  could  be  made  of 
the  burn  distance.  The  samples  burned  from  one  end,  the  cylindrical  surfaces 
being  inhibited  with  a  coating  of  silicone  grease  immediately  prior  to  testing.  Tc 
determine  the  burning  rate  at  a  given  pressure  and  tube  temperature,  three  or 
more  samples  of  each  of  three  different  lengths  were  prepared.  These  were  then 
cemented  to  sample -holders  which  could  be  attached  to  the  end  of  die  push  rod. 

The  burning  times  for  each  of  these  samples  were  then  measured  and  the  average 
burning  rate  computed  from  the  burning  time -sample  length  data.  Figure  2  illus¬ 
trates  the  type  of  photocell  data  that  was  obtained. 

Similar  data  were  also  obtained  with  cylindrical  samples  inhibited  on  the 
end  rather  than  on  the  cylindrical  sides.  These  samples  burned  radially.  The 
diameter  of  these  samples  was  varied,  rather  than  the  length,  and  the  burning 
times  measured  as  a  function  of  radial  burn  distance. 


4 .  Data  Reduction  Method 

An  illustration  of  the  burning  time  data,  plotted  versus  the  sample  length, 
is  shown  in  Figure  3.  As  indicated  by  this  figure,  the  burning  times  can  be  corre¬ 
lated  by  an  equation  of  the  format  -  a  ♦  'bL.  The  reciprocal  of  b,  or  the  slope  of 
the  lines,  represents  the  burning  rate. 

The  best  fit  line  was  determined  by  the  method  of  least  squares  for  each 
set  of  data.  An  estimate  of  the  reliability  of  the  burning  rate  was  then  made  follow¬ 
ing  regression  analysis  procedures  (33).  According  to  these  procedures,  the  esti¬ 
mated  standard  error  in  time  predicted  by  the  best  fit  equation  t  -  a  -  bL  is 

St  -[r  <tj  -  ti)2/(N-2)] 1/2  . (2.1) 

where  tj  is  the  measured  time  and  tj  the  predicted  time  for  a  sample  of  length  Lj, 
and  N  is  the  number  of  data  points.  The  corresponding  standard  error  in  the 
slope  b  is 


Figure  2.  Example  of  photo  cell  data  used  to  compute  burning  rates 
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Figure  3.  Example  of  tube -furnace  data  illustrating  effect  of  thermal 
radiation  on  burning  times. 
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sb  .  V[E  C*i  -  *)2  )1/2  . (2.2) 

1  i 

where 

x  =  Z  xt/N  . (2.3) 

i 


The  95  per  cent  confidence  limits  in  would  then  be  given  approximately  by 

b  =  b  -  2Sb  . (2.4) 


The  corresponding  95  per  cent  confidence  limits  on  r  *  1/S  would  then  be  given 
to  the  same  degree  of  approximation  by 


r  -  ?  1  2Sb  /  (£  )2 


(2.5) 


5.  Experimental  Results 

Table  1  summarizes  the  experimental  burning  rate  data  obtained  during 
this  project.  The  measured  burning  times  and  sample  lengths  for  each  test  are 
tabulated  in  the  Appendix .  Five  different  propellants  were  tested  during  the  course 
of  the  testing  program.  All  of  the  propellants  were  ammonium  perchlorate  oxidized 
composites.  The  propellant  designated  E-107  has  a  polyurethane  binder  with  alum¬ 
inum.  Propellant  UG  consists  of  18  per  cent  polybutadiene -acrylic  acid  binder  with 
82  per  cent  ammonium  perchlorate.  Propellant  AGC  64-1106  has  a  carbaxytermin- 
ated  polybutadiene  binder  with  aluminum .  Propellant  A -13  consists  of  a  poly¬ 
butadiene -acrylonitrile  binder  and  ammonium  perchlorate  in  the  ratio  76/24. 
Propellant  AAP-3318  is  similar  to  AGC  64-1106  but  with  part  of  the  ammonium 
perchlorate  replaced  with  potassium  perchlorate.  The  test  condition  variables, 
in  addition  to  furnace  temperature,  were  sample  configuration  and  furnace  pres¬ 
sure. 


Effect  of  furnace  temperature.  The  effect  of  the  furnace  tube  temperature 
on  the  burning  rates  at  10.3  psia  of  propellants  AGC  64-1106,  UG,  and  E-107  are 
shown  in  Figure  4 .  Approximate  incident  radiant  flux  levels  corresponding  to  the 
measured  tube  temperatures  are  shown  on  this  figure.  These  flux  levels  were 
estimated  using  the  equation 

0  ■  Fl2  «-  '  V4 
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TABLE  1.  Summary  of  Data  Obtained  with  Tube -Furnace 
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Burning  Rate,  inches/second 


Furnace  Temperature  -  °F 
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Figure  4 .  Data  showing  effect  of  thermal  radiation  on  burning  rates 
atl0.3psia. 
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with  F}2  assumed  to  be  0.8,  Ts  assumed  to  be  450°C  (940°F),  and  Tw  the  wall 
temperature.  Recent  attempts  at  measuring  Ts  (5)  for  similar  AP-oxidized  pro¬ 
pellants  have  shown  it  to  be  in  the  range  of  500-740°C  at  pressures  near  100  psia. 
The  temperature  would  be  expected  to  be  lower  at  the  pressures  employed  for 
the  furnace  tests,  assuming  an  Arrhenius  relationship  between  Ts  and  burning  rate. 
Consequently,  the  value  of  450°C  was  selected  as  a  representative  value  in  order 
to  estimate  the  net  radiant  flux .  It  is  noted  that  with  these  assumptions  there  is  a 
loss  of  approximately  0.4  cal/cirr  sec  from  the  burning  surface  when  the  propellant 
is  exposed  to  room -temperature  surroundings. 

The  data  shown  in  Figure  4  indicate  that  the  burning  rates  of  the  UG  and 
AGC  propellants  were  increased  by  approximately  50  per  cent  as  the  furnace  tube 
temperature  was  increased  from  room  temperature  to  1500°F,  while  for  E-107 
propellant  a  50  per  cent  increase  occurred  when  the  temperature  was  raised  from 
room  temperature  to  1000°F.  These  data  clearly  indicate  a  strong  effect  of 
thermal  radiation  on  burning  rate  at  this  pressure. 

Data  for  A -13  and  AAP-3318  propellants  illustrated  in  Figure  5  do  not 
show  the  same  magnitude  effect  as  for  the  other  propellants.  Difficulty  in  igniting 
these  propellants  was  experienced,  however,  and  the  data  are  somewhat  question¬ 
able  . 


Effect  of  sample  configuration.  Initial  tests  with  the  AGC  64-1106  propellant 
•vere  made  with  radial  burning  cylindrical  samples  of  propellant  inhibited  on  the 
ends.  Subsequent  tests  made  with  end -burning  cylindrical  samples,  with  the  sides 
inhibited,  resulted  in  substantially  higher  average  burning  rates  with  the  furnace  at 
room  temperature.  The  room  temperature  mean  rate  for  radial  burning  was  .028 
inches/sec  while  the  end-burning  mean  rate  was  .037  inches/second. 

This  difference  was  thought  to  be  due  to  the  proximity  of  the  cold  stainless 
steel  tube  to  the  burning  surface.  The  minimum  clearance  between  the  burning 
surface  and  the  tube  wall  in  this  configuration  was  0.66  inches,  whereas  in  the  end- 
burning  configuration  the  burning  surface  was  perpendicular  to  the  center  line  of 
the  tube  and  the  motion  of  the  combustion  products  was  unobstructued .  To  test 
this  hypothesis,  tests  were  made  using  the  radial -burning  configuration  with  the 
furnace  tube  replaced  by  a  large -diameter  lucite  tube.  In  these  tests  the  clearance 
between  the  burning  surface  and  the  tube  wall  was  increased  to  2.60  inches.  The 
mean  burning  rate  at  room  temperature  was  observed  to  be  .030  inches/sec  in 
this  configura-tion,  nearly  the  same  as  with  the  small -diameter  tube,  indicating 
that  the  cold  wall  did  not  cause  the  rate  to  be  reduced. 

Figure  6  presents  a  comparison  of  the  data  for  the  different  configurations. 
There  appears  to  be  no  effect  of  the  configuration  when  the  tube  is  heated;  however, 
taking  into  account  the  uncertainty  in  the  data,  the  radial-burning  configuration 
appears  to  result  in  lower  rates  for  this  propellant  when  the  tube  is  unheated.  No 
suitable  explanation  for  this  effect  can  be  given  at  this  time. 
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'I'.sts  with  both  radial  and  end-burning  configurations  were  also  made  with 
the  Ur;  propellant.  For  this  propellant,  the  differences  in  the  mean  rates  were 
within  the  uncertainty  limits,  indicating  no  significant  effect  of  configuration  differ¬ 
ences  . 


Effect  of  pressure.  Figure  7  illustrates  the  effect  of  varying  the  pressure 
in  the  tube -furnace .  This  figure  presents  data  obtained  with  UG  propellant  at 
pressures  of  2.5,  4.9,  and  10.9  psia .  The  room  temperature  (75°F)  data  is  seen 
to  correspond  well  with  conventional  strand  burner  data  for  this  propellant  reported 
by  Lockheed  Propulsion  Company  (5). 

Effect  of  radiation  on  PpL«  A  limited  study  was  made  of  the  effect  of  the 
thermal  radiation  from  the  furnace  tube  on  the  deflagration  limit  of  the  propellant. 
The  results  of  these  experiments  are  also  listed  in  Table  1. 

These  experiments  were  conducted  by  placing  two  valves  between  the  tube 
furnace  and  the  evacuated  tank,  a  metering  valve  and  a  fast -acting  solenoid  valve. 
Prior  to  ignition,  the  pressure  in  the  tube  was  adjusted  to  slightly  less  than  atmos¬ 
pheric  with  the  solenoid  valve  closed.  The  sample  was  then  ignited  and  the  solenoid 
valve  opened,  causing  the  tube  pressure  to  drop  at  a  rate  governed  by  the  netering 
valve.  The  deflagration  limits  were  assumed  to  be  indicated  by  the  photocell  signal 
falling  off.  The  sample  was  withdrawn  from  the  heated  tube  at  this  point,  the  remain 
ing  unbumed  propellant  confirming  that  extinguishment  had  occurred.  To  insure 
that  a  depressurization  effect  was  not  affecting  the  data,  the  experiments  were 
repeated  at  different  metering  valve  settings . 

Figure  8  presents  data  obtained  in  this  manner  for  the  AGC  propellant.  These 
data  indicate  the  Pdl,  was  reduced  approximately  from  4  in.  Hg  to  3  in.  Hg,  absolute 
pressure,  when  the  sample  was  exposed  to  1000°F  tube  walls.  Additional  data  for 
this  propellant  and  UG  propellant  are  listed  in  Table  1. 

Attempts  at  measuring  the  Pql  with  a  tube  temperature  of  1500°F  were 
unsuccessful.  The  samples  always  were  consumed  even  though  the  photocell  indi¬ 
cated  extinguishment  might  have  occurred. 

Conclusions.  These  tests  show  that  quantitative  measurements  can  be  made 
in  the  laboratory  which  show  the  effect  of  incident  thermal  radiation  on  solid  pro¬ 
pellant  burning  rates.  At  sub -atmospheric  pressures,  where  the  burning  rates  are 
low,  radiant  fluxes  corresponding  to  surrounding  walls  at  1500°F  can  cause  the  burn¬ 
ing  rates  to  be  as  much  as  50  per  cent  greater  than  those  measured  under  conven¬ 
tional  strand  burner  conditions .  This  large  effect  should  obviously  be  considered  in 
the  ballistic  design  of  controllable  solid  propellant  motors,  and  the  tube -furnace 
technique  employed  in  this  project  can  be  applied  to  obtain  this  data.  Additional 
work  should  be  done  to  improve  the  accuracy  of  the  data. 
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DATA  COMPILATION  AND  CORRELATION 
1.  Survey  of  Published  Experimental  Data 

This  survey  covers  in  an  historical  fashion  the  work  done  in  the  experimental 
field  of  propellant  extinguishment  via  rapid  depressurization.  Data  from  nine  differ¬ 
ent  sources  are  tabulated  in  this  section.  Table  2  lists  the  nomenclature  used  in 
these  tables. 

References  6,  7,  and  8  describe  the  work  of  Ciepluch  who  performed  the 
original  set  of  experiments  in  this  field.  The  extinguishment  experiments  were 
performed  by  burning  a  slab  of  propellant  in  a  specially  designed  small  motor  that 
contained  two  nozzles .  The  slab  weighed  about  1  pound  and  was  placed  on  the  side 
of  the  motor  so  that  the  gas  flow  was  parallel  to  the  propellant  surface.  The  exper¬ 
imental  technique  consisted  of  first  igniting  the  small  motor  and  permitting  it  to 
reach  a  stable  operating  pressure,  one  nozzle  being  closed.  At  this  time,  the 
second  nozzle  was  opened  by  means  of  an  explosive  bolt,  causing  rapid  depressuri¬ 
zation  .  The  motor  was  then  examined  to  determine  whether  there  was  any  pro¬ 
pellant  left  or  whether  it  had  all  been  consumed.  If  no  propellant  remained,  the 
test  was  classified  as  a  non -extinguishment  test.  If  propellant  remained  in  the 
motor,  extinguishment  was  said  to  have  occurred.  Table  3  summarizes  the  exper¬ 
imental  data  reported  in  these  three  references. 

Marginal  extinguishment  conditions  were  expressed  in  terms  of  a  character¬ 
istic  depressurization  time  above  which  the  grain  burned  out  without  extinguishment. 
This  characteristic  time  was  defined  as  that  necessary  to  depressurize  the  chamber 
to  one -half  its  initial  pressure  and  was  designated  as  tj/2*  Later  investigators  have 
used  both  this  time  and  the  corresponding  average  depressurization  rate  between  the 
initial  pressure,  P*,  and  Pj/2.  Figure  9,  taken  from  reference  7,  illustrates  this 
kind  of  experimental  data. 

The  Ciepluch  technique  was  followed  by  nearly  all  subsequent  investigators 
except  that  in  some  cases  the  propellant  configuration  was  altered  and  different 
methods  were  used  to  increase  the  nozzle  area. 

Reference  9  describes  the  work  conducted  by  Amcel  Propulsion  Company. 

The  extinguishment  portion  of  that  program  was  secondary  to  the  major  objectives 
of  developing  a  controllable  motor,  and  only  one  propellant  was  tested.  Table  4 
lists  the  data  for  the  propellant  tested  in  that  program. 

In  1964,  Aerojet  General  Corporation  began  a  major  program  in  this  area. 
Table  5  summarizes  data  extracted  from  the  reports  describing  the  work  performed 
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Figure  9.  Example  of  method  of  correlating  P  extinguishment 
conditions  (Ref.  7) 
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under  several  different  contracts  (2,  10,  11,  12).  The  test  motor  geometry  used 
for  most  of  the  data  reported  was  very  similar  to  that  used  by  Ciepluch  except  for 
the  fact  that  3"  O.D.,  one-pound,  end-burning  grains  were  used  rather  than  side¬ 
burning  slabs.  Some  data  were  also  obtained  in  large  motors,  however. 

Table  6  summarizes  the  data  obtained  by  Hercules,  Inc  .,  at  their  Bacchus 
Plant  (13).  In  this  investigation  they  studied  the  extinguishability  of  double-base 
solid  propellants  and  used  a  motor  which  was  virtually  identical  to  that  of  Ciep- 
luch's . 


Table  7  and  Reference  14  present  the  data  and  results  obtained  at  the  Univer¬ 
sity  of  Utah .  Their  experimental  procedure  was  unique  in  that  they  used  a  strand  of 
propellant  burning  in  a  very  large  volume  so  that  a  single  nozzle  was  used  to  control 
the  depressurization.  This  nozzle  was  closed  at  the  beginning  of  a  test.  Thus,  in 
this  experimental  technique,  the  pressure  in  the  chamber  would  begin  to  rise  very 
slowly  upon  ignition  of  the  strand.  However,  because  the  strand  burning  area  was 
small  compared  with  the  free  volume,  the  pressure  rise  rate  was  negligible. 

After  stable  combustion  was  realized,  the  nozzle  was  opened,  causing  depressuriza¬ 
tion.  With  this  type  of  experimental  setup,  the  single  nozzle  not  only  controls  the 
depressurization  rate  but  the  final  pressure  if  extinction  does  not  occur. 

The  results  of  a  Stanford  Research  Institute  program  are  described  in 
Reference  15  and  Table  8.  Only  a  limited  amount  of  experimental  data  were 
obtained  in  this  program,  the  effort  being  largely  theoretical.  Although  differing 
in  some  details,  the  experimental  motor  and  technique  were  essentially  that  used 
in  the  Aerojet  program. 

Table  9  and  Table  10,  respectively,  describe  the  results  recently  obtained 
at  the  United  Technology  Center  as  reported  in  References  16  and  17.  In  these 
programs,  a  very  extensive  series  of  propellant  formulations  were  tested  for 
extinguishability  characteristics .  The  extinction  technique  was  similar  to  that 
used  by  Ciepluch.  However,  several  different  grain  configurations  were  used  in 
this  study,  including  strands,  end -burning  grains,  slabs,  and  tubular  internal 
burning  grains.  Because  the  data  from  these  programs  are  presented  more  com¬ 
pletely  than  those  from  other  programs,  Tables  9  and  10  contain  the  data  probably 
of  most  use  to  other  investigators. 

Reference  18  and  Table  11  describe  the  results  of  an  experimental  program 
carried  out  at  Brigham  Young  University.  The  technique  used  was  essentially  that 
of  Aerojet  except  the  end -burning  grains  were  only  1-1/4  inches  in  diameter. 

Because  it  has  become  apparent  that  the  data  as  reported  in  Reference  18  were 
not  sufficient  for  many  purposes,  the  data  contained  in  the  present  survey  have 
been  expanded  over  that  originally  presented. 


20 


Table  12  describes  the  results  obtained  in  the  program  of  Atlantic  Research 
Corp.  (19).  The  experimental  technique  used  was  essentially  that  used  at  the  Uni¬ 
versity  of  Utah  in  that  an  extremely  large  volume  of  gas  made  the  use  of  a  primary 
nozzle  unnecessary  and  a  single  nozzle  controlled  the  entire  blowdown.  References 
20-24  also  contain  additional  experimental  data  describing  the  extinguishment  of 
burning  solid  propellants .  However,  either  the  data  were  taken  in  a  much  different 
manner  than  described  above  or  are  not  reported  in  sufficient  detail  to  permit  com¬ 
parison  with  the  data  listed  herein. 
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TABLE  2.  Nomenclature  Used  in  Survey  Tables  (Tables  3-12) 

The  area  of  the  secondary  nozzle  divided  by  the  primary  nozzle  area 

P/po 

The  exponent  in  the  burning  rate  law  r*apn  at  PQ 

The  initial  depressurization  rate  when  the  secondary  nozzle  opens 

The  ambient  pressure 

Chamber  pressure  prior  to  the  opening  of  the  secondary  nozzle 
The  deflagration  limit  of  the  propellant 
Burning  rate  of  the  propellant  at  Pj 

The  time  required  for  the  chamber  to  depressurize  to  Pj  /2 
Adiabatic  flame  temperature 


TABLE  3.  Extinguishment  Data,  NASA  Lewis  Research  Center 

Investigator:  Carl  C .  Ciepluch 
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TABLE  4.  Extinguishment  Data,  Amcel  Propulsion  Co. 


Propellant 

PBO- 

Extinguishment 

P,  10^  psi/sec 

38 

Pi,  psia 

1000 

rj,  in/ sec 

.27 

"i 

.65 

Pa,  psia 

14.7 

Reference:  Second  Annual  Report,  RRL-PT-4-64-52,  September,  1964. 


TABLE  5.  Extinguishment  Data,  Aerojet  General  Corporation 
Investigator:  R.  Lou  and  L.  Landers 
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TABLE  6.  Extinguishment  data,  Hercules,  Inc. 


Propellant 

CYH 

E]C 

VHX 

VHY 

VHW 

VHV 

VCP 

VIR 

(Modified  Double  Base) 

Pdi,  psia 

1.45 

8.2 

2.3 

1.9 

1.9 

4.7 

2.3 

3.0 

Pi ,  psia 

512 

366 

420 

360 

380 

340 

400 

400 

380 

r| ,  in/sec 

.434 

.378 

.41 

ni 

.55 

.55 

Pa,  psia 

12.5 

0.5 

0.5- 

12.5 

.5 

.5 

.5 

Extinguish. 

P ,  10^  psi/sec 

60 

10 

9.5 

9.5 

12 

15 

7 

7 

7 

Reference:  Final  Report,  Contract  No.  AF04(694)-127  WS-133A,  1965  PSP  Task  8 
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TABLE  7.  Extinguishment  Data,  University  of  Utah 


Propellant 

G 

AH 

GB 

UA 

AF 

F 

UF 

Binder 

PBAA 

PBAA 

PBAA 

PBAA 

PBAA 

PBAA 

PU 

Cone .  % 

18 

25 

17 

15 

15 

18 

20 

Oxidizer 

Ap 

Ap 

Ap 

Ap 

Ap 

Ap 

Ap 

Cone  .  % 

82 

75 

80 

73 

73 

80 

80 

Size, 

200/15 

15 

200/15 

15 

200/15 

200/15 

200/15 

Additive 

Carbon 

Black 

CU0202* 

cuo2o2 

cuo2o2 

cuo2o2 

Cone .  % 

2 

2 

2 

2 

Tf »  °K 

2625 

1988 

2459 

2598 

Pdl,  psia 

.7 

.92 

.49 

2.2 

Pif  psia 

90-175 

70-150 

90-175 

90-175. 

90-175 

90-175 

rj ,  in/sec 

.16-. 22 

.17-. 23 

.16-. 22 

.37-. 49 

i 

.28-. 42 

.15-. 26 

ni 

.5 

.56 

.50 

.55 

.38 

.78-. 31 

Pa,  psia 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

Exting. 

dlnP/dt 

50 

50 

50 

225 

600 

700 

11 

Reference:  Final  Report,  Contract  AFOSR  67-1901,  September,  1966. 

*CU0202  denotes  copper -chromite  manufactured  by  Harshaw  Chemical  Co. 


TABLE  8.  Extinguishment  data,  Stanford  Research  Institute 


Propellant 

PU174 

PU193 

PU185 

Binder 

PU 

PU 

PU 

Cone .  % 

17.5 

20 

18.5 

Oxidizer 

AP 

AP 

AP 

Cone.  % 

80 

80 

80 

Size.yu. 

Additive 

Ethyl 

Siloxane 

Fe2°3 

Cone.  % 

2.5 

1.5 

Tf,  °K 

2600 

Pdl,  psia 

3.3 

1.0 

3 

Pi ,  psia 

500 

470 

750 

ri  ,  in/sec 

.19 

.32 

.37 

ni 

.6 

.5 

.5 

Pa,  psia 

7.8 

5.2 

6.0 

Extinguishment 

P,  103  psi/sec 

19 

14 

27 

Reference:  QTR  No.  5,  Contract  NAS  7-389,  February,  1967. 


TABLE  9.  Extinguishment  Data,  United  Technology  Center 

Investigator;  C.  E.  Jensen 
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Reference:  Final  Technical  Report,  Contract  NAS  1-6601,  UTC  2243-FR,  November,  1967. 
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Investigator:  G.  D.  Jensen 
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TABLE  11.  Extinguishment  data,  Brigham  Young  University 


Propellant 

A -13 

A-14 

A-15 

A  -16 

A-17 

A -18 

Binder 

PBAN 

PBAN 

PBAN 

PBAN 

PBAN 

PBAN 

Cone .  % 

24 

24 

24 

24 

24 

24 

Oxidizer 

AP 

AP 

AP 

AP 

AP 

AP 

Cone .  % 

76 

76 

76 

76 

76 

76 

Size,/< 

80 

15 

80 

15 

80 

15 

Additive 

-- 

CUO^Ofc  CLOp2 

LiF 

LiF 

Cone .  % 

-- 

1 

1 

1 

1 

Tf,  °K 

2100 

2100 

2100 

2100 

2100 

2100 

Pdi,  psia 

.40 

.56 

.20 

.45 

.10 

.80 

Pj  ,  psia 

91 

164 

400 

153 

280 

285 

320 

169 

157 

Tj  ,  in/sec 

.105 

.143 

.23 

.24 

.39 

.31 

.67 

.121 

.181 

ni 

.51 

.51 

.51 

.61 

.61 

.53 

.56 

.42 

.57 

Pa,  psia 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

Extinguishment 

1 1/2 ,  msec  11 

10 

7.5 

7.2 

7.2 

6.4 

1.9 

7.3 

6.6 

Anf/^nj 

1 .99 

2.92 

4.50 

2.39 

2.78 

4.40 

5.15 

4.99 

2.57 

Reference:  AIAA  Journal,  6,  292-297  (1968) 


TABLE  12.  Extinguishment  data,  Atlantic  Research  Corp. 

Investigators:  G.  Von  Elbe,  E.  McHale 


Prop . 1  lant 

Arcite 

PBAA 

PBAA-A1 

Binder 

PVC 

PBAA 

PBAA 

Cone.  % 

20 

30 

28.5 

Oxidizer 

AP 

AP 

AP 

Cone.  % 

80 

70 

66.5 

Size, 

80 

80 

Additive 

A1 

Cone.  % 

13.4 

Pi,  psia 

100 

200 

200 

in/sec 

.174 

.109 

.178 

ni 

• 

c/I 

.50 

.42 

Pa.  psia 

14.7 

14.7 

14.7 

Extinguishment 
tl/2’  msec* 

20 

51 

13.4 

♦Computed  from  data  presented 


Reference:  AIAA  Journal,  6,  July  1968,  1417-1419. 
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2.  Correlations  of  Experimental  Data 

A  first  attempt  at  correlating  the  data  discussed  in  the  previous  section  util 
ized  the  following  equation,  which  was  derived  independently  by  Von  Elbe  (27)  and 
Paul,  et  al ,  (28): 


or 


1  r^ 

-  (dlnP/dt)ext  .  T  —  (3.1) 

•  <«  V'2”  ■  77" 


where  (dlnP/dt)ext  is  the  critical  logarithmic  depressurization  rate  for  extinguish¬ 
ment,  r  is  the  burning  rate,  n  is  the  burning  rate  exponent,  *  is  the  thermal  diffus- 
ivity,  and  A  is  an  empirical  correction  factor.  Figure  10  presents  the  results  of 
this  correlation.  To  obtain  this  plot  it  was  assumed  that  cc  =  0.00025  in^/sec  for 
all  propellants.  This  assumption  was  necessary  since  data  for  «  are  not  available 
for  most  of  the  propellants.  It  was  further  assumed  that  the  pressure  decayed  expon¬ 
entially  in  all  cases.  The  theory  leading  to  Equation  (3.1)  predicts  that  extinguish¬ 
ment  occurs  during  the  pressure  decay  transient  when  this  equation  is  satisfied.  Thus 
the  value  of  r  employed  in  the  correlation  should  be  that  corresponding  to  steady -state 
at  the  pressure  when  dlnP/dt  satisfied  Equation  (3.1).  Since  only  average  initial 
decay  rates  are  reported,  the  time  during  decay  when  extinguishment  occurred  is 
not  known.  If  the  decay  is  exponential,  the  value  of  dlnP/dt  is  constant  during  the 
depressurization  and  the  average  initial  value  will  apply  approximately  throughout 
the  decay.  The  minimum  value  of  the  right-hand  side  of  (3.1)  would  be  that  corre¬ 
sponding  to  the  final  pressure,  since  P  and  hence  r  would  then  have  their  minimum 
values.  For  this  reason,  the  correlation  was  attempted  (Figure  10)  using  the  burning 
rates  corresponding  to  the  final  pressure.  This  attempted  correlation  was  clearly 
unsatisfactory. 

Figure  11  presents  the  results  of  a  second  correlation  attempt.  In  this  figure 
the  depressurization  half-time,  which  is  proportional  to  dlnP/dt  for  an  exponential 
decay,  is  plotted  versus  the  initial  steady  burning  rate .  Although  the  data  show  a 
trend  consistent  with  Equation  (3.1),  the  scatter  is  such  that  this  correlation  was  also 
not  satisfactory. 

As  a  result  of  the  study  discussed  in  Section  IV,  it  became  apparent  that  the 
motor  L*  is  an  important  variable  that  must  be  accounted  for  in  any  successful 
correlation  of  extinguishment  data.  The  importance  of  this  variable  is  illustrated 
by  the  data  shown  in  Figure  12,  which  were  taken  from  Reference  10.  These  data 
show  that  die  marginal  initial  depressurization  rate  for  extinguishment  for  a  given 
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Figure  10.  Results  of  attempt  to  correlate  marginal  P 
extinguishment  with  burning  rate  exponent 


Time  to  50%  -  t>  -  Milliseconds 


2 


Marginal 

Initial 

Depressurization 

Rate 

(103)psi/sec 


□  Swing  Nozzle  Motor 
15"  Grain,  L*  *  70" 

□  Swing  Nozzle  Motor 
S"  Grain,  l  *  100" 


□  Window  Motor 
L*  %  2000" 


Pressure,  Psia 


1000  121 


Figure  12. --Data  taken  from  Reference  17  showing  the  effect 
of  L*  on  marginal  1*  extinguishment. 
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propellant  can  vary  by  about  a  factor  of  5  depending  on  the  motor  used  for  testing. 

The  critical  rate  appears  to  increase  with  decreasing  initial  L*  of  the  test  motor. 

As  discussed  in  Section  IV,  the  marginal  extinguishment  conditions  predicted 
theoretically  were  nicely  correlated  without  regard  for  depressurization  rate  using 
the  parameters  L^*  and  r^.  These  parameters  are  the  characteristic  length  of  the 
motor  and  the  burning  rate  that  would  be  calculated  assuming  steady-state  ballistics 
apply  using  the  fully  opened  nozzle  area,  or  using  Kn  *  Ajj/Anf .  Figure  13  presents 
the  results  of  the  correlation  attempt  which  was  made  following  this  approach . 

Unfortunately,  most  of  the  literature  cited  previously  in  this  section  does  not 
contain  sufficient  information  to  permit  application  of  this  method  of  correlation.  A 
notable  exception  is  the  recent  data  reported  by  the  United  Technology  Center  (16,  17), 
and  most  of  the  data  points  shown  in  Figure  13  were  taken  from  this  reference.  In 
addition  to  this  data,  one  point  was  extracted  from  results  reported  by  Stanford 
Research  Institute  (15),  one  point  from  full-scale  motor  testing  carried  out  by  Aero¬ 
jet  General  Corp.  (25),  and  one  from  previous  tests  at  BYU.  Also  included  in  tills 
figure  are  theoretical  points  and  lines  illustrating  the  predicted  effect  of  flame  tem¬ 
perature.  These  theoretical  results  are  discussed  in  Section  IV. 

It  is  important  to  note  that  the  correlations  shown  in  Figure  11  and  13  both 
suggest  that  the  lower  the  burning  rate  of  a  propellant,  the  more  easily  it  is  extin¬ 
guished. 
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IV 


PARAMETRIC  DESIGN  STUDIES 


1 .  Objectives  and  Approach 

The  ballistics  engineer  must  make  one  basic  calculation  in  the  design  of  a 
single -chamber  controllable  motor  having  stop-restart  capability.  This  is  the  cal¬ 
culation  of  the  minimum  area  to  which  the  nozzle  must  be  opened  in  order  to  cause 
the  propellant  to  extinguish.  Three  criteria  have  been  proposed  to  calculate  the 
required  area  increase:  (1)  the  depressurization  rate  resulting  from  the  area  change 
be  greater  than  some  critical  value  (P  criterion),  (2)  the  predicted  steady-state 
chamber  pressure  resulting  from  the  change  be  less  than  some  critical  value  (Pql 
criterion),  and  (3)  the  L*  resulting  from  the  change  be  less  than  some  critical  value 
(L*  criterion). 

If  there  exists  a  unique  deflagration  limit,  P^*  or  minimum  pressure  below 
which  a  propellant  will  not  burn,  the  nozzle -area  increase  necessary  to  reach  this 
PpL  can  be  computed  in  a  straightforward  manner  using  conventional  steady -state 
ballistics.  There  appears  to  be  a  strong  effect  of  the  thermal  environment  on  the 
PDL  (26)>  however,  and  proven  methods  have  not  yet  been  developed  for  measuring 
in  the  laboratory  PDL's  which  apply  inside  full-scale  rocket  motors.  In  addition  to 
this  problem,  the  required  area  changes  computed  solely  on  the  basis  of  the  Pql 
criterion,  without  regard  to  motor  L*  effects  or  P  effects,  tend  to  be  larger  than 
those  resulting  from  employing  these  alternative  criteria.  Thus,  the  area  change 
computed  solely  on  the  basis  of  a  limiting  Pql  may  dictate  a  larger  and  heavier  nozzle 
than  is  actually  required. 

The  objective  of  this  phase  of  the  program  was  to  carry  out  parametric  calcu¬ 
lations  for  the  minimum  nozzle -area  change  that  will  lead  to  extinguishment  taking 
both  f*  and  L*  effects  into  account.  The  results  of  these  calculations  were  then  to  be 
graphically  presented  in  a  form  that  would  be  useful  for  preliminary  design  calcula¬ 
tions  . 

The  approach  that  was  followed  was  to  develop  a  mathematical  model  of  the 
combustion  chamber,  coupling  the  gas  dynamic  equations  to  the  thermal  dynamic 
equations  describing  the  combustion  process,  and  then  to  compute  predicted  depres¬ 
surization  transients,  varying  the  operating  parameters . 

2.  Analytical  Model 

Conventional  ballistics  equations  for  a  solid  propellant  rocket  motor  apply 
strictly  only  to  the  steady-state  since  steady-state  burning  rates  are  assumed  to 
apply.  The  conditions  leading  to  extinguishment  are  very  non-steady  and  since 
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experimental  data  for  non-steady  burning  rates  are  difficult  to  obtain  and  not  readily 
available,  ballisticians  have  employed  theoretical  models  of  the  combustion  process 
to  predict  non -steady  behavior.  In  particular,  the  rather  simple  models  of  Von 
Elbe  (27),  Paul  (28),  and  Cohen  (34)  have  been  used  in  controllable  motor  develop¬ 
ment  programs  (25,  26). 

These  models  lead  to  the  P  criterion  mentioned  in  Section  III, 


-  (dlnP/dt)ext  £  r2/  A  n  «. 


(4.1) 


and,  as  shown  by  Cohen  (34),  to  the  L*  extinguishment  criterion: 


(L*)ext  *  nf2C*«/r2 


(4.2) 


As  indicated  in  Section  HI  of  this  report,  Equation  (4.1)  has  limited  useful¬ 
ness  in  correlating  experimental  data  for  different  propellants .  This  is  most  likely 
due  to  the  fact  that  the  theoretical  basis  of  these  equations  employs  rather  restrictive 
assumptions .  In  particular,  the  rate  of  heat  conduction  into  the  solid  is  assumed  to 
occur  instantaneously,  resulting  in  a  temperature  profile  and  heat  storage  that  is  a 
function  only  of  the  instantaneous  value  of  the  pressure.  They  also  assume  the  sur¬ 
face  temperature  to  be  a  constant . 

More  recent  theories  for  extinguishment  have  been  proposed  by  Horton  (18), 
Wooldridge  and  Marxman  (15),  and  Summerfield,  et  al.  (35).  In  these  theories  non¬ 
steady  heat  conduction  is  accounted  for  using  numerical  analysis  techniques.  The 
recent  theories  are  similar  in  many  respects  to  theories  originally  developed  for 
application  to  combustion  instability,  for  example,  the  Denison-Baum  theory  (29). 

Predicted  responses  of  the  burning  rate  to  imposed  pressure  transients  is 
illustrated  by  the  curves  in  Figure  14.  This  figure  shows  the  type  of  response 
predicted  by  the  Von  Elbe  equation  to  a  ramp  or  step  change  in  pressure  relative 
to  the  type  of  response  predicted  by  the  Denison-Baum  combustion  theory.  Using 
the  Von  Elbe  equation,  the  assumption  of  instantaneous  heat  transfer  results  in  the 
prediction  of  a  discontinuous  burning  rate  response.  The  lag  and  overshoot  response 
predicted  by  the  Denison-Baum  theory  is  certain  to  be  more  nearly  the  kind  of 
response  that  actually  occurs . 
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Figure  14.  Comparison  of  burning  rate  response  predicted  by 
Yon  Fibe  equation  with  Denison -Baum  theory. 
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a .  Combustion  Model 

The  combustion  model  employed  in  this  analysis  was  developed  at 
BYU  under  an  AFOSR  Grant  (30).  It  relies  on  a  previous  theory  developed 
by  Denison  and  Baum  (29);  however,  their  theory  was  modified  to  predict 
transient  deviations  from  empirical  steady-state  burning  rates.  The  basic 
features  of  this  model  are  described  in  the  following. 

The  burning  process  is  assumed  to  be  represented  by  the  one -dimen¬ 
sional  temperature  profile  shown  in  Figure  15.  The  solid  is  divided  into 
finite -difference  elements  as  is  also  illustrated  in  this  figure,  and  the  follow¬ 
ing  energy  balance  taken  on  an  element  of  the  solid  is  assumed  to  describe 
the  transient  heat  conduction. 


dTi  . 

■  ■  -  ot 

'  Ti-i-2T|  -Tt-i 

♦  r 

Ti.rTi-i 

dt 

• 

2Ax 

The  instantaneous  burning  rate  is  coupled  to  the  surface  temperature, 
pressure,  and  flame  temperature  by  the  following  equations: 

r  .  ;p  exp  [<-Es/R>  <1/TS  -  1/TSiP>]  (4.4) 

r  .  rp  (T/T,)”*1  exp  [-(^/R)  (1/T(-1/T(  p)]  (4.S) 

and  the  heat  flux  supplied  to  the  solid  from  the  combustion  process  is  given 
by 

dT  _ 

<  -k  aTW  *  r rs  [Cscr,-T.) .  Cg  <Tf-Tf)]  (4.6) 

In  these  equations 

rp  =  the  steady -state  burning  rate  at  pressure  p 
Es  «  Activation  energy  for  surface  reaction 
R  *  Universal  gas  constant 
_  Ts  ■  Instantaneous  surface  temperature 
TSn  Steady  state  surface  temperature  at  p 
_  Tf  =  Instantaneous  flame  temperature 
Tf  p  =  Steady -state  flame  temperature  atp 
’n  -  dlnydlnPat  P 

Eg  -  Activation  energy  for  gas -phase  reaction 
( s  »  Density  of  solid 
Cs  «  heat  capacity  of  solid 
Cg  -  heat  capacity  of  gas 
Too  "  conditioning  temperature  of  propellant 


Geometrical  Model  of  Burning  Propellant 


Figure  15.  Diagram  showing  temperature  profile  during  combustion  and 
method  of  separating  solid  into  finite  elements. 


so 


The  surface  flux  is  incorporated  into  the  finite -difference  energy 
balance  equations  by  the  following  equation  defining  T0,  the  temperature 
at  a  fictitious  station  immediately  above  the  surface. 


T0  -T2  ^2Ax(-k^-)XBO+  (4.7) 


Special  procedures  for  reducing  the  number  of  finite  difference  ele¬ 
ments  needed  for  accurate  calculations  are  discussed  in  Reference  30. 

b .  Interior  Ballistic  Model 

The  model  for  the  solid  propellant  combustion  process  has  been  com¬ 
bined  with  equations  describing  the  transient  ballistics  of  a  solid  rocket 
motor.  Equations  expressing  conservation  of  both  mass  and  energy  inside 
the  motor  cavity  have  been  derived  previously  (31).  In  the  present  study  these 
equations  have  been  used  in  the  following  form: 


*  *  v  (PA c)  <V?C1)  [<VTc>  -  wn|  (4.8) 

^  *  <TC/P)  &  -  <Tc/rc>  <P,/P)  <TC/TC)  ((r/F,)  -  W„|  (4.9) 

where 

rc  -  l\/r  2C* 


y*l 


► 

2 

i  *  * 

Y-l 

lV41 

wn  *  <TC/TC)1/2  (p/p)  (VS) 


Prediction  of  the  instantaneous  pressure  versus  time  requires  the 
simultaneous  numerican  integration  of  these  equations  along  with  those  pre¬ 
sented  in  the  previous  section  for  the  solid  burning  rate.  This  is  accom¬ 
plished  with  an  integration  scheme  based  on  the  Adams  method  (32).  The 
Fortran  coding  for  carrying  out  these  calculations  using  the  BYU  Librascope 
L-30S5  computer  is  presented  along  with  an  example  of  the  computed  output 
in  the  Appendix . 
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3  .  Comparison  of  Predicted  and  Experimental  Transients 

The  theoretical  model  which  has  been  described  above  is  based  on  several 
simplifying  assumptions  which  can  be  justified  only  if  the  model  is  useful  for 
correlating  experimental  data  and  predicting  transient  behavior  as  motor  and 
operating  parameters  are  varied.  Evaluating  the  model  for  these  purposes,  of 
course,  depends  on  the  availability  of  suitable  data.  Although  the  developed  model 
allows  for  the  effect  of  external  radiation,  this  effect  was  not  included  in  the  calcu¬ 
lations  performed  here. 

The  model  requires  as  input  data  all  of  the  usual  parameters  necessary  to 
calculate  the  steady-state  operating  pressure  of  a  rocket  motor  plus  several  addi¬ 
tional  motor  and  propellant  parameters.  The  additional  motor  parameters  are 
the  initial  L*.  the  change  in  venting  area  that  takes  place  to  terminate  thrust, 
and  the  time  required  to  change  the  venting  area.  The  additional  propellant 
parameters  are  the  surface  temperature  at  some  specified  burning  rate,  Tsr,  the 
steady-state  flame  temperature,  the  gas  activation  energy,  Eg,  the  solid  phase 
activation  energy,  Es,  the  ratio  of  the  solid  and  gas  phase  heat  capacities,  and  the 
thermal  diffusivity  of  the  solid.  All  of  these  additional  parameters  can  be  readily 
obtained  except  Tsr,  Es,  and  Eg.  These  become  the  curve  fitting  or  correlating 
parameters.  The  theory  is  most  useful  only  if  one  set  of  these  parameters  for  a 
given  propellant  can  be  used  to  predict  transient  behavior  as  operating  conditions 
are  varied. 

The  method  that  has  been  used  to  select  these  parameters  is  to  assume  a 
Tsr  typical  of  those  measured  with  very  fine  thermocouples  (600°C)  and  then  to 
adjust  Es  and  Eg  as  required  to  fit  existing  experimental  data.  For  convenience, 

Es  and  Eg  have  been  taken  to  be  equal,  so  that  only  a  single  parameter  is  adjusted 
to  reach  agreement  with  a  given  set  of  data .  This  approach  has  been  satisfactory 
for  correlating  two  sets  of  existing  data  where  the  initial  pressure  and  ambient 
pressure  were  varied;  i.e.,  one  set  of  constants  for  a  given  propellant  was  suffi¬ 
cient  to  predict  the  effects  on  extinguishment  conditions  observed  experimentally 
as  these  operating  parameters  were  varied. 

The  experimental  data  referred  to  were  taken  from  a  series  reported  by 
Horton  (18)  and  from  a  different  series  reported  by  Wooldridge  (15)  of  the  Stanford 
Research  Institute.  Horton  determined  the  critical  depressurization  rate  necessary 
to  cause  extinguishment  of  a  propellant  designated  as  A -13  at  initial  pressures  of 
170  and  310  psia.  The  ambient  pressure  for  these  tests  was  12 .5  psia.  The  SRI 
researchers  determined  the  critical  vent  area -nozzle  area  ratio  necessary  to 
extinguish  a  propellant  designated  as  PU-269  at  ambient  pressures  of  15  psia  and 
165  psia,  die  initial  pressure  being  held  constant  at  approximately  500  psi .  Table 
13  summarizes  the  propellant  input  data,  including  the  activation  energies  finally 
selected  to  correlate  the  data,  and  Table  14  summarizes  die  motor  input  data 
and  lists  the  observed  and  predicted  results.  As  mentioned  above,  Tsr  was 
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TABLE  13.  Summary  of  Propellant  Characteristics 
Employed  in  Extinguishment  Predictions 


Propellant  Designation 

A-13 

PU-269 

Binder 

PB AN/  Epoxy  24% 

Polyurethane  20% 

Oxidizer 

Am.  Perch.  76% 

Am.  Perch.  80% 

Input  Propellant  Parameters 

Burning  Rate,  100  psi,  in. /sec. 

0.33 

.33 

Exponent,  1000  psi 

0.41 

.40 

Flame  Temperature,  °K 

2500 

2910 

Surface  Temperature,  1000  psi,  °C 

600 

600 

Density,  lbs/in'* 

.061 

.062 

Cp/Cyi  Gas 

1.22 

1.22 

C*,  ft/sec 

3070 

4500 

<Cp)  solid/  <C«)  gas 

Diffusivity,  iir/sec 

0.75 

0.75 

.00025 

.00025 

Conditioning  Temperature,  °C 

25 

25 

E8,  Kcal/mole 

20 

21 

Eg,  Kcal/mole 

20 

21 

TABLE  14.  Comparison  of  Predicted  and  Experimental  Extinguishment  Data 
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assumed  to  be  600°C  for  both  propellants .  The  magnitude  of  E  ■  Eg  ■  Es  was 
selected  by  a  trial  and  error  approach,  noting  that  increasing  E  causes  extinguish¬ 
ment  to  occur  more  readily.  The  activation  energies  finally  selected  were  those 
that  correctly  predicted  marginal  extinguishment  at  both  extremes  of  the  operating 
conditions,  i.e.,  chamber  pressures  of  170  and  310  psia  in  die  case  of  A-13  pro¬ 
pellant,  and  ambient  pressures  of  15  and  165  psia  in  the  case  of  PU-269  propellant. 

Figures  16  and  17  illustrate  predicted  transients  for  four  of  the  cases  sum¬ 
marized  in  Table  14.  Figure  16,  along  with  the  data  in  Table  14,  shows  that 
increasing  the  vent  area/nozzle  throat  area  ratio  from  3 .9  to  5.6  in  the  SRI  exper¬ 
iments  causes  a  decrease  in  the  depressurization  half-time  from  8  to  4  milli¬ 
seconds,  causing  the  burning  rate  to  drop  to  levels  where  extinguishment  would  be 
expected  to  occur  (burning  rates  less  than  about  .005  inches/sec  are  not  observed 
to  occur  during  steady -state  and  would  not  be  expected  during  transients).  Thus  a 
blow -down  half-time  tj/2  °*  between  4  and  8  milliseconds  is  predicted  as  marginal 
for  extinguishment  for  this  propellant  at  this  particular  set  of  operating  conditions. 
This  prediction  is  in  reasonably  good  agreement  with  the  experimental  critical  tj/2 
between  3  and  5  milliseconds  observed  by  SRI  under  conditions  near  those  assumed 
for  the  prediction . 

A  reasonably  severe  test  of  the  theory  is  to  determine  whether  or  not  it  can 
be  employed  to  predict  the  observed  effects  of  increasing  the  ambient  pressure  from 
15  psia  to  165  psia  without  an  adjustment  being  required  in  the  effective  activation 
energies,  Es  and  Eg.  After  several  trials  in  which  various  values  of  Es  and  Eg 
were  used  as  input  data,  a  value  of  21  kcal/mole  was  found  to  satisfy  the  observed 
effect  with  reasonable  accuracy.  Using  this  value,  a  critical  tj/2  of  4-8  milli¬ 
seconds  was  predicted  for  an  ambient  pressure  of  15  psia,  as  mentioned  above,  and 
a  critical  tj/2  °f  0*6  -1.2  milliseconds  was  found  for  an  ambient  pressure  of  165 
psia;  the  latter  value  compares  well  with  the  observed  value  of  0.6  -  1.0  milliseconds. 

Figure  17  shows  two  of  the  predicted  transients  made  in  this  series  of  calcu¬ 
lations  which  illustrate  the  strong  effect  of  ambient  pressure.  With  ambient  pressure 
of  165  psia  extinguishment  falls  to  occur  even  though  t^2  was  near  2  milliseconds. 
With  ambient  pressure  of  15  psia,  however,  extinguishment  occurs  with  tj/2  near 
8  milliseconds .  Thus  the  marginal  tjy2  or  dP/dt  is  predicted  to  depend  strongly  on 
the  ambient  pressure . 

A  different,  though  somewhat  less  severe,  test  of  the  model  is  to  determine 
if  the  observed  effect  of  varying  the  initial  pressure  can  be  predicted  with  a  single 
set  of  activation  energies.  This  test  was  made  for  die  BYU  experiments  of  Table  14. 
As  listed  in  Table  14,  the  theory  predicts  the  critical  tj/2  is  increased  from  8-18 
milliseconds  to  12-24  milliseconds  if  the  initial  pressure  is  reduced  from  309  to 
167  psia.  These  figures  bracket  the  critical  times  observed  in  the  tests  carried 
out  at  BYU.  This  agreement  was  achieved  assuming  Es  -  Eg  ■  20  kcal/mole,  as 
listed  in  Table  13 . 


Time  -  milliseconds 


Figure  16.  Predicted  termination  transients  for  input  parameters 
corresponding  to  PU-269  propellant  (Ref.  15). 


4.  Results  of  Parametric  Calculations 


Because  of  the  large  number  of  input  variables  that  are  contained  in  the 
mathematical  model  described  above,  the  effects  of  varying  each  parameter  were 
not  exhaustively  studied  during  this  program.  The  calculations  were  limited  to 
those  showing  the  depressurization  rates  that  occur  under  certain  limiting  condi¬ 
tions,  the  overall  pressure  drop  that  results  if  extinguishment  does  not  occur, 
and  the  effect  on  marginal  extinguishment  conditions  of  varying  initial  pressure, 
initial  L*,  and  nozzle  opening  time. 

a .  Limiting  Depressurization  Rates  and  Overall  Pressure  Drop 

For  the  special  case  of  an  instantaneous  opening  of  the  throat  area, 
Equation  4.8  reduces  to  Equation  4.11  for  the  initial  depressurization 
rate. 


1  dP  12vr2C*  An  .  x 

Pdt  Lj  [a^ 


(4.11) 


where  C*  has  the  units  of  feet/second  and  Lf  is  the  initial  L*  in  inches. 
Figure  18  shows  die  initial  fractional  depressurization  rate  versus  the 
instantaneous  fractional  increase  in  nozzle  throat  area,  with  the  initial 
L*  as  a  parameter. 


If  the  nozzle  is  opened  slowly  enough,  steady -state  ballistics  apply 
and  the  depressurization  rate  is  given  by 


(4.12) 


for  a  propellant  whose  steady -state  burning  rate  can  be  approximated  by 
r  =  a  I  .  Figure  19  shows  a  plot  of  these  limiting  fractional  depressuriza¬ 
tion  rates  versus  the  fractional  rate  of  nozzle  area  increase. 


Figures  18  and  19  are  useful  for  bracketing  the  depressurization 
rates  that  may  be  computed  more  exactly  using  the  complete  set  of  equa¬ 
tions  of  Section  3 . 


A  third  set  of  curves  representing  limiting  conditions  are  presentee 
in  Figure  20.  The  curves  in  this  figure  represent  solutions  to  the  equation 


(4.13) 


which  relates  the  fractional  change  in  steady-state  pressure  to  the  fractional 


Figure  19.  Predicted  depressurization  rates  assuming  the  burning  rates  and  ma* 
flow  in  and  out  of  the  chamber  are  at  equilibrium. 
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change  in  nozzle  throat  area  for  a  propellant  following  an  r«aPn  burning  rate 
law  and  for  which  C*  is  independent  of  pressure. 

Figures  18  and  20  can  be  employed  to  point  out  an  important  depres¬ 
surization  relationship,  the  relationship  between  initial  depressurization 
rate,  initial  L*  and  overall  pressure  drop.  As  an  illustration  of  this  rela¬ 
tionship,  consider  the  following  experiment.  Five  different  motors  are 
fired  having  initial  L*’s  of  100",  500",  1000",  1500",  and  2500".  They  are 
at  the  same  initial  pressure,  500  psia,  and  their  nozzles  are  instantaneously 
opened  to  the  point  where  the  initial  depressurization  rate  is  20,000  psi/ 
second.  If  extinguishment  does  not  occur,  the  new  steady-state  pressures, 
according  to  Figures  18  and  20,  will  range  from  345  psia  for  the  100"  L* 
motor  to  9  psia  for  the  2500"  L*  motor,  as  shown  in  the  following  table. 


Initial  L* 
inches 

Fractional  Nozzle 
Increase  for 

(dP/dt)i  ■  20,000  psi/sec 

Final  Steady -state 
Pressure 
n  =  0.6 

100 

0.16 

345  psia 

500 

0.80 

115 

1000 

1.60 

46 

1500 

2.40 

24 

2500 

4.00 

9 

These  final  steady-state  pressures  of  course  apply  only  for  a  Vielle  pro¬ 
pellant  (?  =  aP0),  but  the  same  trend  shown  by  this  exercise  should  occur 
for  a  real  propellant;  i.e.,  the  nozzle -area  changes  that  will  result  in  the 
same  initial  dP/dt  can  result  in  drastic  differences  in  the  overall  pressure 
drop  depending  on  differences  in  the  initial  L*. 

If  extinguishment  results  not  only  as  a  result  of  achieving  a  critical 
dP/dt  but  also  because  of  a  large  overall  pressure  drop,  these  calculations 
show  that  differences  in  motor  size  can  have  an  important  effect  on  marginal 
extinguishment  conditions.  The  calculations  discussed  in  the  following  sec¬ 
tion  indicate  that,  from  the  ballistic  designer’s  point  of  view,  a  critical  over¬ 
all  pressure  drop  is  the  most  important  criterion  for  achieving  extinguish¬ 
ment.  They  show  further  that  expressing  extinguishability  of  a  propellant  in 
terms  of  a  critical  dP/dt  alone  is  rather  meaningless . 

b.  Basis  for  Transient  Model  Calculations 

Parametric  calculations  showing  the  effect  on  marginal  extinguishment 
conditions  of  varying  initial  pressure,  initial  L*,  and  nozzle  opening  time 
were  all  made  using  a  single  set  of  model  propellant  parameters.  The  param¬ 
eters  for  the  model  propellant,  designated  propellant  Y,  are  listed  in  Table  15. 
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TABLE  15.  Properties  of  Theoretical  Propellant  Y 


Burning  Rate  @  1000  psia,  in/sec  0.400 

Burning  Rate  Exponent  @  1000  psia  0 .400 

Flame  temperature,  °K  3000 

Burning  Surface  Temperature  at  1000  psia,  °C  600 

Density,  (  ,  lbm/in3  .060 

Ratio  of  Specific  Heats,  r  1 .20 

Characteristic  Velocity,  C*,  ft/sec  3000 

Surface  Activation  Energy,  Es,  Kcal/mole  20 

Gas  Activation  Energy,  Eg,  Kcal/mole  20 

Solid  Heat  Capacity/Gas  Heat  Capacity,  C8/C  0.50 

© 

Solid  Thermal  Diffusivity,  « ,  in^/sec 


00025 
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Its  burning  rate  is  shown,  plotted  versus  pressure,  in  Figure  21.  The 
corresponding  Kn  curve  is  also  shown.  Since  the  burning  rate  exponent 
(dlnr/dlnP)  for  most  propellants  is  not  constant  over  wide  variations  in 
pressure,  a  Summerfield  burning  rate  law  (35)  was  assumed. 

The  reference  steady-state  motor  operating  conditions  were  assumed 
to  be  500  psia,  which,  with  the  assumed  propellant  properties,  corresponds 
to  a  motor  Kj,  of  186.3.  The  ambient  pressure  in  all  cases  was  assumed  to 
be  1  psia . 

To  study  the  effect  of  initial  pressure  variation,  calculations  were 
made  at  100  psia  and  1000  psia  in  addition  to  the  reference  pressure  of  500 
psia.  The  initial  L*  was  varied  between  10  inches  and  5000  inches,  and  the 
nozzle  opening  time  was  varied  between  0.2  milliseconds  and  413  milli¬ 
seconds.  For  each  combination  of  LJ,  and  nozzle  opening  time,  the 
minimum  nozzle  area  change  resulting  in  extinguishment  was  calculated. 

The  results  of  all  of  these  calculations  are  summarized  in  Table  16. 

c .  Effect  of  Varying  Initial  Pressure 

The  results  of  the  calculations  made  to  illustrate  the  effect  of  vary¬ 
ing  initial  pressure  are  presented  in  Figure  22,  following  the  procedure 
which  has  been  employed  by  others  (6)  to  correlate  P  extinguishment  data. 
For  these  calculations,  the  initial  L*  was  1000  inches  and  the  opening  time 
was  0.2  milliseconds.  The  marginal  initial  depressurization  rate  is  shown 
to  increase  uniformly  with  increasing  initial  pressure,  just  as  observed 
experimentally,  provided  a  single  motor  and  propellant  configuration  is 
used  for  each  test . 

d.  Effect  of  Varying  Initial  L* 

Figure  23  presents  the  results  of  calculations  in  which  the  initial 
L*  was  varied  between  100"  and  5000",  the  nozzle  opening  time  being  0.2 
milliseconds  and  the  initial  pressure  being  500  psia  in  all  cases .  The 
data  are  plotted  in  a  similar  manner  to  those  in  Figure  22  except  that 
interpolated  marginal  initial  depressuration  rate  is  plotted  instead  of  the 
points  for  individual  calculations.  The  marginal  depressurization  rate  is 
shown  to  increase  uniformly  with  decreasing  LJ.  The  similarity  between 
these  results  and  the  experimental  data  of  UTC  shown  in  Figure  12  is  quite 
significant. 

e.  Effect  of  Varying  Nozzle  Opening  Time 

Increasing  the  nozzle  opening  time,  the  overall  nozzle-area  change 
and  other  parameters  remaining  constant,  naturally  resulted  in  a  reduction 


TABLE  16.  Summary  of  Parametric  Study  Calculations,  Propellant  Y 
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Figure  22 .  Predicted  effect  of  initial  pressure  on  p  for 
marginal  extinguishment. 
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Figure  23.  Predicted  effect  of  initial  L*  on  P  for  marginal  extin 
guishment  conditions. 
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in  the  depressurization  rates.  The  surprising  result  of  these  calculations, 
however,  was  that  the  marginal  overall  nozzle -area  change  was  affected 
only  slightly.  Figure  24  presents  a  comparison  of  area -change  versus 
initial  L*  for  opening  times  of  0.2  milliseconds  and  20  milliseconds,  the 
initial  pressure  being  500  psia.  This  figure  shows  that  the  marginal  nozzle 
area  change  is  only  slightly  affected,  if  at  all,  by  die  change  in  opening 
time. 


Figure  25  illustrates  the  predicted  effect  of  the  opening  time  on  the 
depressurization  rate.  The  depressurization  rates,  expressed  as  dlnP/dt, 
differ  by  approximately  a  factor  of  four  until  the  pressure  reaches  100  psia. 

At  this  pressure,  which  corresponds  to  a  time  just  prior  to  the  nozzle  being 
completely  open,  the  dlnP/dt  for  the  slowly -opened  case  suddenly  increases 
to  near  that  of  the  rapidly -opened  case.  This  rapid  change  is  apparently  due 
to  transient  combustion  effects .  The  important  point  illustrated  by  this  figure 
is  that  the  initial  depressurization  rate,  and  also  the  depressurization  rate 
that  would  be  predicted  using  steady -state  burning  rates,  is  of  little  value  to 
the  designer  in  predicting  marginal  extinguishment.  For  both  cases  illus¬ 
trated,  the  nozzle-area  ratio  is  near  the  marginal  ratio  below  which  extin¬ 
guishment  does  not  occur. 

f.  Correlation  of  Marginal  Extinguishment  Conditions 

After  several  attempts  using  different  approaches,  it  was  discovered 
that  all  of  the  results  of  this  parametric  study  could  be  correlated  in  terms 
of  the  motor  L*  following  opening  of  the  nozzle,  LJ,  and  the  burning  rate 
that  would  result  from  the  nozzle-area  change,  ff,  provided  steady-state 
ballistics  apply.  Thus  Lj>  «  V.Anf  and  for  a  propellant  with  a  burning  rate 
law  r  =  aP11,  rf  -  rj  (Anj/Anf)n^  .  Figure  26  presents  all  of  the  data  from 
Table  16,  taken  without  regard  for  nozzle -opening  times,  plotted  as  L| 
versus  r{.  The  line  Lf  *  1.3/fj  ,  which  has  the  same  form  as  Equation 
(4.2),  is  shown  to  provide  a  reasonable  correlation  of  the  marginal  extinguish¬ 
ment  conditions.  Thd«,  it  appears  that  if  the  mathematical  model  employed 
in  these  computations  accurately  describes  die  real  transient  process 
occurring  in  a  rocket  motor,  the  ballistic  designer  can  predict  marginal 
extinguishment  nozzle  sizes  without  regard  to  depressurization  rates. 

The  correlation  shown  in  Figure  13  of  Section  III  presents  additional 
support  both  for  the  use  of  the  Denison -Baum  theory  and  the  correlative 
approach  of  Figure  26  for  predicting  extinguishment  conditions.  As  was 
previously  noted,  Figure  13  contains  both  theoretically  predicted  marginal 
extinguishment  conditions  and  the  motor  data  available  with  sufficient  sup¬ 
porting  data.  Using  the  propellant  properties  of  Tables  13  and  15  for  the 
respective  flame  temperatures,  the  marginal  conditions  predicted  by  the 
Denison -Baum  approach  are  shown  in  Figure  13.  As  delineated  by  the 
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Figure  25.  Predicted  depressurization  curves  showing 
the  effect  of  nozzle -opening  times. 
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Final  Burning  Rate,  inches/second 

Figure  26.  Correlation  of  parametric  study  results  in  terms  of 
steady-state  L*  and  burning  rate  following  the  nozzle 
area  change. 
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dashed  lines,  the  theoretical  and  experimental  data  show  a  reasonable 
separation  based  on  flame  temperature.  Hence  the  theoretical  and  exper¬ 
imental  data  show  good  agreement .  Based  on  the  limited  data  for  compar¬ 
ison,  the  correlative  approach  of  Figures  13  and  26  appears  sound. 
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RECOMMENDED  DESIGN  PROCEDURES 

As  a  result  of  the  studies  carried  out  during  this  project,  certain  recom¬ 
mendations  can  be  made  for  improving  the  procedures  which  have  been  employed 
to  predict  marginal  extinguishment  conditions  in  single -chamber  controllable  motors. 
These  recommendations  fall  into  two  categories:  (1)  those  having  application  to  pre¬ 
liminary  design  calculations,  and  (2)  those  having  application  to  experimental  devel¬ 
opment  work  leading  to  final  designs. 

It  is  cautioned  that  these  recommendations  are  of  a  tentative  nature.  Fur¬ 
ther  parametric  studies  should  be  made  to  investigate  the  importance  of  additional 
operating  variables  such  as  conditioning  temperature,  irregularities  in  the  steady- 
state  burning  rates,  and  the  thermodynamic  and  kinetic  properties.  More  impor¬ 
tantly,  additional  experimental  data  should  be  obtained  to  more  fully  establish  the 
reliability  of  the  combustion  model  that  has  been  employed. 

It  is  further  pointed  out  that  these  recommendations  involve  only  the  predic¬ 
tion  of  marginal  extinguishment  conditions,  and  that  the  problem  of  reignition  follow¬ 
ing  extinguishment  has  not  been  considered. 

1 .  Preliminary  Design 

For  carrying  out  preliminary  design  studies  that  include  determining  opti¬ 
mum  operating  conditions  based  on  weights  of  inert  motor  components,  it  is  recom¬ 
mended  that  the  practice  of  applying  the  Von  Elbe -Paul  theory  to  define  a  critical 
depressurization  rate  be  discontinued.  The  simple  criterion  resulting  from  the 
study  described  in  Section  IV,  and  applied  to  existing  data  in  Section  III,  namely, 
that  extinguishment  will  occur  when 


Lf  s  B  (5.1) 

is  recommended.  It  has  shown  promise  of  being  applicable  regardless  of  the  rate 
at  which  the  nozzle  area  is  changed,  the  size  and  configuration  of  the  motor,  and 
the  initial  steady-state  operating  conditions  that  are  imposed.  Furthermore,  it  is 
extremely  easy  to  apply,  requiring  only  standard  steady-state  ballistic  calculations. 

The  magnitude  of  the  constant  B,  which  may  be  interpreted  as  a  parameter 
characteristic  of  a  given  propellant,  will  typically  range  from  0.2  to  1.3  when  LJ 
is  expressed  in  inches  and  ff  in  inches  per  second,  as  is  shown  by  the  data  of  Fig¬ 
ure  13. 
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Even  though  B  for  different  propellants  may  range  from  0.2  to  1.3,  the 
sensitivity  of  the  predicted  marginal  nozzle  -area  change  to  variation  in  this  param¬ 
eter  is  rather  small,  as  illustrated  in  Figure  27.  In  this  Figure  the  magnitude  of 
the  product  Lf  r2  is  plotted  versus  the  corresponding  nozzle -area  change  ratio. 

The  initial  pressure  for  Aerojet  General  Corp.  formulation  AAB  3220  is  assumed 
to  be  520  psia  and  the  initial  L*  is  taken  at  9.55  inches.  For  this  case,  the  pre¬ 
dicted  marginal  nozzle -area  ratio  ranges  only  from  2.8  to  4.0  as  B  ranges  from 
0.22  to  1.3.  Thus,  even  assuming  the  full  range  uncertainty  in  B,  ignoring  flame 
temperature  effects,  the  uncertainty  in  the  marginal  nozzle  area  ratio  is  approxi¬ 
mately  t20%. 

2 .  Experimental  Development  Work 

The  analytical  study  carried  out  during  this  project  showed  that  character¬ 
izing  the  extinguishability  of  a  propellant  in  terms  of  a  critical  initial  average 
depressurization  rate  provides  design  information  of  very  limited  usefulness. 

This  kind  of  marginal  extinguishment  criterion  can  be  directly  applied  to 
full-scale  motor  design  only  if  the  initial  L*  and  pressure  of  this  motor  are  the 
same  as  those  used  in  obtaining  the  experimental  data .  A  criterion  showing  promise 
of  much  greater  generality  is  given  by  Equation  (5.1).  Therefore,  it  is  recommended 
that  the  experimental  test  program  be  planned  with  the  objective  of  determining  the 
magnitude  of  the  constant  B  in  this  equation . 

Using  this  approach,  testing  would  be  carried  out  in  much  the  same  manner 
as  used  before  to  obtain  P  data .  However,  the  important  data  to  be  recorded  are 
the  burning  areas,  nozzle  areas,  combustion  chamber  volumes,  and  steady-state 
pressures . 

Since  the  motor  L*  appears  in  the  transient  ballistic  equations  along  with 
C*,  in  the  group  L*/  2C*,  attention  should  be  given  to  the  effective  C*  that  results 
from  small  motor  tests .  If  the  C*  is  lower  in  these  firings  than  anticipated  in  the 
full-scale  motor,  then  the  B  used  for  the  full-scale  design  should  be  increased,  or 


bFS  "  Bss  <C*FS/C*SS> 


where  the  subscript  FS  denotes  full  scale  and  the  subscript  SS  denotes  sub-scale. 

This  approach  emphasizes  the  importance  of  obtaining  accurate  steady- 
state  burning  rate  data,  particularly  at  low  pressures.  A  technique  to  obtain  labor¬ 
atory  data  more  representative  of  that  expected  in  the  full-scale  motor  is  therefore 
recommended.  The  tube -furnace  approach  described  in  detail  in  Section  n  provides 
useful  data  and  is  a  recommended  technique. 
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Initial  Throat  Area/Final  Throat  Area 

.O 

Figure  27.  Example  showing  the  variation  of  the  product  L*r  with 
nozzle -area  change  illustrating  effect  of  uncertainty  in 
marginal  extinguishment  criterion. 
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The  Denison -Baum  combustion  theory  which  has  been  employed  in  this 
study  assumes  a  rather  simplified  geometrical  model  to  represent  the  actual 
combustion  process.  However,  on  the  basis  of  the  limited  comparisons  of  pre¬ 
dicted  and  transient  behavior  presented  herein,  it  shows  promise  as  an  effective 
curve -fitting  model.  It  has  been  shown  to  correctly  predict  the  effects  of  varia¬ 
tions  in  the  motor  operating  parameters .  It  can  likely  be  used  as  well  to  correlate 
the  effects  of  variations  in  the  propellant  parameters  that  can  be  obtained  directly, 
such  as  flame  temperature,  steady-state  burning  rate,  and  burning  rate  exponent. 
Furthermore,  it  provides  a  basis  for  describing  the  transient  burning  properties 
of  the  propellant  in  terms  of  effective  activation  energies .  It  is  therefore  recom¬ 
mended  that  an  effort  be  made  to  characterize  the  transient  combustion  behavior 
of  propellants  in  terms  of  this  model. 
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NOMENCLATURE 

1 .  Experimental  Data 

An  Nozzle  throat  area 

B  Constant 

C  Heat  capacity 

C*  Characteristic  velocity 

E  Activation  energy 

Burning  area /throat  area 
k  thermal  conductivity 

L*  Characteristic  length 

n  Burning  rate  exponent  •  dlni/dlnP 

P  Pressure 

P  Rate  of  depressurization 

Prtl  Deflagration  limit 

r  Burning  rate 

T  Temperature 

t  time 

ll/2  Time  to  depressurize  to  Pj/2 

X  distance 

Thermal  diffusivity 
Y  Ratio  of  specific  heats 
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Empirical  constant 
Density 

(  "  )  Bar  denotes  steady -state 


a 

c 

f 

g 

i 


Flame,  final  value 
Gas 


Initial  value 


p  Value  at  pressure  p 

s  Surface,  solid 

2 .  Computer  Program 

LENGTH  Sample  length  or  radius 

TIME  Measured  bum  time 


ERROR  Difference  between  measured  time  ana  time  predicted  by  best 
fit  equation 

RBAR  Least -squares  best  fit  burning  rate 

ERROR  95%  confidence  limit  for  RBAR  (Equation  2.5) 
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APPENDIXES 


APPENDIX  I 

TABULATION  OF  BURNING  RATE  DATA  OBTAINED 
IN  THE  TUBE -FURNACE  APPARATUS 


PRECEDING  PAGE  BLANK 


F-107 


£-107 


TEMP  21  IN  HC 


TlNt 

30.10. 
28 . 20  J 

27.10. ’ 
27.20c 
17 , 90ki 
16.80. 
17.50*. 
16.40. 
16.00, 
17.200 
12.40L 
18.00. 
13 . 80k. 
13.40c 
12.30c, 
12.50  c’ 


ERROR 

2.417 

0.817 

-0.502 

-0.402 

0.039 

-1.060 

-0,360 

-1.460 

-1.860 

-0.660 

-0.201 

2.398 

0.898 

0.790 

-0.301 

•0.101 


RAOUL  RM 
LENGTH 
0.3950 
0.3950 
0.3950 
0.3960 
C.24tf0 
0.2400 
0.2400 
C .2400 
0.2400 
0.2400 
0.1570 
0.1570 
0.1570 
0.1570 
0.1570 
0.1570 


RADIAL 
LENGTH 
0.3950 
0.3950 
0.3950 
0.3950 
0.24pe 
0.2400 
0.2400 
0.2400 
C.2400 
6.2400 
C .  2400 
0.1570 
0.1570 
0.1570 
0.1576 
0.1570 
0.1570 
0.1570 
0.1570 


RBAH  s 
FRRUR  s 


i:00  OEG  F 
TIME 
16.600 
17.500 

17.300 

18.300 
11 . 20c 

10.700 

10.70J 

10.70k 

10.70c 

10.20k 

10.300 
6.60J 

7.300 
7.60c 
8*600 
6.80c 
7.900 
7  •  60c 
8.50v 

9BAR  * 
ERROR  s 


0.01578 

0,60166 


21  IN  HG 

error 

-0.6/4 

0.225 

0.625 

1.025 

0.309 

-0.190 

-0.190 

-0.190 

-0.190 

-0.690 

-0.590 

-0.871 

•0.171 

0.126 

1.128 

-0.671 

0.428 

0.126 

1.028 

0.02427 

0.00166 


PRECEDING  PAGE  BLANK 


TEMP  5 

In  hg 

7/16/6 V 

LENGTH 

Tint 

ERROR 

e .3950 

22.00. 

1.037 

0.395k, 

19.20.. 

-1.76' 

3.2400 

1.4.20. 

-0.00-5 

0.2400 

15.10 

0.894 

3.2400 

15.80. 

1 .594 

0.2400 

13.00. 

-0 . 40-3 

3.1570 

9.30.; 

-1 .  26« 

0.1570 

11,50... 

0.91 

0.1570 

9.70., 

—  0  »  H  6  b 

3.1570 

10.50. 

•0.06b 

RRAR  s 

0.02294 

FRROR  = 

0.00448 

ibZZ  DEG 

F  10 

IN  HG  RAD 

LENGTH 

TIME 

error 

0.3950 

12.20. 

—  0 . 306 

3.3950 

12.50* 

-0.006 

0.3950 

12.60* 

0.293 

0.2400 

8.20C 

0.052 

0.2400 

8.00c 

-0.147 

‘ 0.2400 

8.30* 

0.152 

0.1570 

5.60J 

-0.212 

0.1570 

5.60* 

-0.212 

0.1570 

6.200 

0.367 

RRAR  s 

0.03555 

ERROR  a 

0.00221 

uu 


IJG 


hm 


TFMP  1 

IM  HG 

RADIAL 

UG 

RM  TEMP  21 

IN  HG 

RADIAL 

length 

TlMf. 

ERROR 

LENGTH 

TIME 

error 

Z . 3950 

15.50, 

-0.044 

0.3950 

9.20* 

•0.044 

3.3950 

15.80. 

0.255 

2.3950 

9.25. 

0.005 

£.3950 

14.70. 

-0 . 844 

0.3950 

9.25. 

0.005 

3.2400 

9.70, 

-0.195 

0.2400 

5.60* 

-0.066 

3.2400 

10.80  > 

0.904 

0.2400 

5. 75;. 

0.081 

3.24pt0 

11.00. 

1.104 

0.2400 

5.75. 

0.061 

0.1.570 

6.9e. 

0.029 

0.1570 

3.901 

0.146 

0.1570 

5.80. 

-1.070 

0.1570 

3.70* 

-0.053 

0.1570 

6.10. 

-0.770 

0.1570 

3.60.’ 

-0.153 

0.15/0 

7.50. 

0.629 

RRAR  s 

0.04334 

RRAR  * 

0.02743 

ERROR  s 

0.00125 

FrtRUR  s 

0.00360 

• 

UG 

1002  DEC 

E  21 

IN  HG  KADI 

LENGTH 

TIME 

ERROR 

0.3950 

8.70* 

0.076 

133*  n£G 

f  1 ? 

IN  HG  RADIAL 

0 • 3950 

8.60J 

-0.023 

LENGTH 

TlMf 

ERROR 

0.3950 

8.50. 

.0.123 

2.395b 

12.00. 

-0.710 

0.2400 

5.50* 

0.03  3 

0.3950 

11.80 

-0.910 

0.2400 

5.50* 

0.033 

3.3950 

l4.0i\ 

1.269 

0.24?e 

5.60  *T 

~0.133 

0.2400 

9.00* 

0.351 

0.1570 

3.80.. 

0.023 

0,2400 

8.70, 

0.051 

0.1570 

3.70*’ 

-0.076 

0.2400 

9.20, 

0.551 

0,1570 

3.70. 

-0,076 

0.1570 

6.50. 

0.026 

0.1570 

5.90, 

-0.573 

RBAR  s 

0.04910 

0.1570 

6.40 

-0.073 

ERROR  «" 

0.00142 

RRAH  s  0.03815 
CRROR  s  0.00/21 


106:  n£& 

F  21 

IN  HG 

UG 

1586  DEC  F 

21  IN 

HG  RADI 

ENGTH 

TIME 

ERROR 

LENGTH 

TINE 

ERROR 

0. 29*50 

7 . 50  j 

0.156 

0.3950 

7.50 

0.404 

0.2750 

7.2P^ 

0.265 

0.3950 

6.70. 

-0.399 

2.2780 

7.00 

0.065 

2,3950 

7.16. 

0.004 

3.1650 

3.20. 

-0,276 

0.2400 

4.4P. 

-  0 , 3  7  6 

0.1000 

3.30, 

-0,074 

"0.2400 

4796 

f'.Tfi 

0.1000 

3 ,5P 

0.125 

0.2400 

5.00 

0.221 

6.1806 

4.96. 

-0.101 

0.1570 

3.70. 

0.161 

0.19P0 

5.00. 

-0.205 

0.1570 

3.40. 

-0.138 

0.1930 

5.40. 

0.133 

0.4650 

9.86. 

0.221 

RRAH  s 

0.06600 

0.4680 

9.80. 

0.16* 

FRROR  r 

0.01011 

0.49?0 

11.30, 

-0.048 

0.5150 

11.60. 

-0.216 

0.^900 

11.13. 

- 0 . 208 

UG 

Rm  tehp  21 

IN  HG 

STRAND 

RRAR  s 

0.04915 

LENGTH 

TIME 

ERROR 

FRROR  s 

0,00168 

2.1060' 

- 4720. 

0.26? 

0.2980 

4.00, 

0.23* 

0.1080 

3.50. 

-0,480 

0.2770 

7.50. 

-0.060 

DEG  F 

21  IN  HG 

strand 

0.2700 

7.60. 

0.179 

ength 

TIME 

ERROR 

0.2760 

7.30., 

-0.247 

0.T970 

3.20. 

0.064 

0.4810 

11 .70, 

-0.2(90 

0.T970 

3.06. 

-0,135 

0.4300 

11.20w 

0.382 

S.:970 

2.90  j 

-0.235 

6.4650 

11.50, 

-0.000 

0.1606 

4 . 56, > 

0.194 

0.1600 

4.30 

-tr.045 

RBAR  s 

0.04700 

0.1600 

4.40, 

0.054 

FRROR  = 

0.00307 

C . 3900 

9.00. 

0.23/ 

0.3900 

9.00u 

0.23/ 

0.39P0 

9.00. 

0.237 

0.4800 

10.70, 

0.209 

0.4800 

10.40. 

-0.090 

0.4800 

9.80. 

•0.690 

PRAR  s  0.05200 
FRROR  =  0,00280 


GC  RH 


AGC 


AGC 


TEMP  X0  I K 

HG  RADIAL  AGC 

1.00  DEG  F 

2l  IN  HG 

RADIAL 

L  iN'JTH 

T  !  H£ 

ERROR 

LENGTH 

TIHE 

ERROR 

0.3950 

22. 0C. 

-0.0iO 

3.3950 

10.20c 

-0.461 

e.395/ 

20.50. 

-1.539 

0.3950 

10.702 

0 . 03o 

C.3950 

23.53. 

1.460 

0.3950 

11.10c 

0 . 43o 

0.3950 

22.30. 

0.26/ 

0.2400 

6.202' 

-0.24; 

0.2400 

13.50. 

-0.901 

0.2400 

6.60c 

0.152 

0.2400 

14.70. 

0.298 

0.2400 

6.50c 

0.052 

C .  24p0 

14.60. 

0.198 

0.1570 

4.10. 

-0.09/ 

0.1570 

10.66. 

0.288 

3.1570 

4.30c 

0.109 

C .  1570 

10.10. 

-0.211 

0.1570 

4.20. 

0.009 

3.1570 

10.50. 

0.186 

RRAH  X 

0.0367/ 

RBAN  x 

0.02029 

FRROR  s 

0,0024/ 

error  * 

6.00216 

AGC 

150*.  DEG  F 

21  IN 

HG  RADIAL 

16:0  LEG 

F  1.. 

IN  HG  RAD 

LlNGTH 

TIME 

ERROR 

LENGTH 

T  1 ME 

ERROR 

0.3950 

10.102 

0.561 

0.3950 

17.00c 

•0.626 

0.3950 

9.70  2 

0.161 

0.3950 

18.00. 

0.373 

0.3950 

8.90c 

— 0 , 638 

3.3956 

17.70. 

0.073 

0.2400 

6.10c 

" -0.014 

3.2406 

11.60. 

-0.093 

0.2400 

6 . 00c 

-0.114 

C . 2400 

12.30. 

0.60o 

0.2400 

■  '  6".nrj 

- .0,114 

0  •  1570 

8.10. 

-0.416 

0.1570 

4 . 40  j 

0.119 

3.1570 

8.60. 

0.08  5 

0.1570 

4.40Li 

0.119 

0.1570 

4.20^ 

-0.08/ 

PRAR  x 

0.02612 

FRROR  x 

0.00236 

RBAR  r 

0.04526 

errQr  I 

0. 00471 

1500  DFG  F 

1. 

IN  Hi.  RADI  At 

-  •  . . 

LENGTH 

TIHE 

ERROR 

0.3950 

17.20. 

0.434 

0.3950 

15.50. 

-1.265 

0.3950 

16.90. 

0.134 

3.24c*0 

11.50. 

0.698 

3.24P0 

10. 9J. 

/.29rt 

> 

0.2400 

1 1 . 40„ 

0.793 

0.1570 

7.00. 

-0.36/ 

3.1570 

6.60. 

-0. 700 

3 .1570 

7. 20.. 

-0.100 

0.1576 

7.10. 

-0.20/ 

PBAH  * 

0.02514 

FRROR  X 

0.00264 

AGO 

RM  TEMP 

21  IN  Hr, 

SlRANIi 

AGC 

1500  DEG  F 

2t  IN  HG 

STRAND 

LENGTH 

TIME 

ERROR 

LENGTH 

TIME 

ERROR 

0.6000 

18.70, 

-0.110 

0.0970 

3.20, 

-0,06-> 

• 

0.6000 

10.90. 

0.064 

0.0970 

3.10c 

-0.1ttb 

0.6000 

18.50. 

-0.310 

0 . 1600 

4.80. 

0,11a 

r 

0.4000 

13.20, 

-0.212 

Z ,1600 

4.90c 

0.214 

C  •  4070 

14 . 10t 

0.667 

0,1600 

5.10. 

0.414 

G  •  4000 

13.60. 

0.187 

0.3900 

8.90- 

-0.B94 

0 . 2000 

7.70 

-0.31b 

0.3900 

9.20. 

-0 ,5VA 

O.2CP0 

0.00. 

-0.01b 

0.4800 

12.00. 

0.205 

0.4800 

12.10-j 

0.30b 

RBAR  = 

0.03705 

0.4800 

12.30c 

0 . 50b 

ERROR  = 

0.00222 

RHAR  = 

0.04501 

error  = 

0.00363 

AGC 

10-C  deg  f 

21  IN 

HG  STRAND 

LENGTH 

1 1  ME 

ERROR 

0.5000 

14. 9k. 

-0.066 

AGC 

RM  TEHP  21 

IN  HG  PLASTIC  TUBE 

G.5000 

14.40, 

-0.566 

LENGTH 

TIME 

ERROR 

0.5070 

14.70. 

-0.266 

0.3950 

12.80, 

0.109 

0.3070 

10.40. 

0.700 

0.3950 

13.00, 

0.309 

G.3070 

10.20. 

0.500 

0.3950 

12.60c 

-0,090 

0.3000 

10.30, 

0.600 

e.2500 

7.60fc 

-0 . 38* 

0.1000 

4.20. 

-0.233 

07250%' 

'  "7.60. 

—  0 . 3tte 

- 

2.1070 

4.2P, 

-0.233 

0.2500 

8.20, 

0.211 

0.1000 

4.00. 

-0,433 

0.2500 

7.90, 

-0.08* 

0.2500 

7.80, 

-0.18H 

RH  AR  s 

0.03797 

0.1570 

5,00v 

0.02/ 

error  * 

0.00297 

0.1570 

5.10- 

0.127 

0.1570 

5.10  c 

0.12/ 

0.1570 

5.20. 

0.22  7 

RRAR  * 

0 . 03083 

ERROR  s 

0 . 0014  6 

AGC 

RM  TfHP 

21  IN  HG 

strand 

ELASTIC  TUBE 

LENGTH 

TIME 

ERROR 

0.1370 

5.80. 

-0.369 

0.1 3*0 

6.40 

0.230 

0.1370 

6.00, 

-0.169 

3,1370 

6.20 

0.030 

C . 2130 

9.30 

0.366 

3.2130 

8.6C, 

-0,313 

3.2136 

8.80. 

-0.113 

0.2130 

9.00, 

0.086 

G .  2870 

10 . 6  <6  v 

-0.525 

3.2870 

11.50* 

0.371 

0.2800 

11.50. 

0.371 

0.2800 

11.30, 

0.171 

0.5400 

19.80 

0.075 

C  ,5400 

19.50. 

-0.224 

Rfi  Art 
ERROR 


0.03024 

0.0011b 


A -15  S7W.NL  1  0t  * 
LINO  IH 
2.0  5ry; 

:.b5ru 

e. 3«ez 

:.34?e 

0.34*0 

2.19*0 

2.19*50 

2.39Kb 


A-13  STRAND  15 
LENGTH 
0.55*0 
0.55*0 
0.5500 
0.34*0 
2.3400 
C.34*0 
0.34ee 
3.19*50 
S  •  19*0 
C.l9*0 


A-13  STRAND  RH 
L‘ NGTH 
0.55P0 
0.5500 
0.5500 
C . 3400 
0.340b 
0 . 34*0 
0.19*0 
C • 1950 
C.1950 


nf  G  F 

21  ]N  HO 

T  I  Hf 

ERROR 

17.30 

-0.043 

16.70., 

-0.643 

17.1;).. 

-0.243 

12.40. 

0.5*3 

12.30 

0.4*3 

13.00.. 

1.193 

7.60. 

-0.303 

7.50. 

-0.403 

7.50. 

-0.403 

PBAR  s 

0.03792 

r RROR  s 

0.00438 

2  DEG  F 

21  IN 

TlMt 

error 

15.10. 

-0.729 

15.4gw 

-0.42* 

15.33. 

-0.529 

11 . 80 j 

1.104 

11.90. 

1.204 

11.70 

1.804 

11.20c 

0.504 

6.30. 

.0.715 

6.30. 

-0.715 

6.P0. 

- 1 .015 

PRAR  x 

0.04027 

ERROR  x 

0.00723 

TEMP 

21  IN  HO 

TIME 

error 

18.10. 

0.11/ 

17.90.. 

•0.002 

17.30. 

•0.602 

12.90. 

0.262 

12.50. 

-0.137 

14.10. 

1.462 

8.80c 

-0.146 

8.40. 

-0.546 

8.70. 

-0.246 

PRAR  s 

0.0392m 

ERROR  s 

0.00469 

AAP-3318  STRANnS 
I  LNGTH 
8.1 9*5(6 
0.J95& 
0.1950 
3.3980 
0. 3900 
0.3980 
0.5860 
0.5880 
0.5880 


AAP-3318  STRANDS 
LENGTH 
0.1950 
0.1950 
0.1950 
C  .3900 
0.3900 
0.3980 
0.5580 
0.5580 
0.5580 


Hh  Tt HP 

25  IN  Ht, 

TIME 

ERRO* 

10.400 

-0.000 

9. 40. 

-x.002 

9.20* 

-1.200 

18.00.. 

2.363 

16.60. 

0.96/ 

16,70.. 

1.06/ 

19.50, 

-1.45/ 

21.50. 

0.540 

19.70. 

-1.257 

R8AH  = 

0.03723 

ERROR  s 

0.0081o 

1000  DEG 

F  25  IN 

TIME 

ERROR 

6.80. 

0.126 

7.30. 

0.628 

6.9  C: 

0,228 

12.00, 

•0.6// 

12 . 50(o 

•0.1// 

11. 4e, 

-1.2/7 

17 . 30k. 

-0.551 

18.00, 

0.146 

19 . 40w 

1.546 

RBAR  3 

0.03246 

ERROR  3 

0.00412 

AAP-3318 


RADI  A| 

length 
s.2«pe 
8. £400 
2.?4p0 
0.3950 
C . 3950 


1530  DEG  T 
TIME 
6.80. 
6.9^ 
7.20. 
10.800 
10.20. 


25  IN  hg 
ERROR 
-0.165 
-0.066 
0.233 
0.300 
-0.2yv 


RBAR  a  0.043do 
ERROR  3  0,006/5 


APPENDIX  II 


FORTRAN  CODING  AND  EXAMPLE  OF  COMPUTER  OUTPUT 
OF  BALLISTICS  MODEL  CALCULATION 


PRECEDING  PAGE  BLANK 


PRECEDING  PAGE  BLANK 


COMPUTER  PROGRAM  FLOW  SHEET 


019.T-10  COATFS 


FORTRAN  DECK  *ma!N 


t 


10/24/69 


C 

C 


1 

1 

1 

1 

1 

1 

1 

2 

3 

4 

5 

10 

6 

7 

13 

8 

9 

14 

10 

11 

12 

43 

13 

14 

15 

16 

17 

18 

44 

19 

20 

45 

21 

22 

23 

40 

24 

25 

26 

30 

27 

41 

27 

• 

28 

42 

28 

• 

28 

• 

28 

• 

28 

• 

28 

• 

20 

• 

20 

• 

20 

• 

29 

30 

80 

30 

• 

30 

• 

31 

32 

33 

34 

34 

35 

35 

36 

17 

36 

58 

50 

19 

MOTOR  thrust  TERMINATION 
R.L. COATES  REVISED  OCTOBER  1969 

COMMON  TPRNT(4)» Jl, J2, J3,NPTS» I. 

to; pa, ab.xkn,anr,topen, elstar.eps, 

R8AR,RMIN,XNBaR,TFBAR,TSBAR.RH0S»GAM,CSTAR. 

ES. EG, CSCG.FEXT, ALPHA, 

OAPGAM,PCR!T,GAMFl,GAMF2,GAMF3f TCHAM, TSOI » BB» DC , TC I . 

X. P, R.OLPDT, PATE, PRESS. TF.TC.TS.TL AST, DEPTH, AN. XNP, 
TARLP(A) .TABLR(8),RP 
COMMON  J,TIME(100> .PRES (10&). TEMP (100), BRA TE< 100) 

DIMENSION  NAME ( 20  ) 

REAO(10)0 

J«0 

READ  (5.13)NAME 
FORMAT  (20A4) 

WRITE(6,16)NamE 
FORMAT  (1H1.20A4) 

N»12 

READ (5, 43)  TPRNT(1),TPRNT(2), Jl, J2.J3.NPTS 
FORMAT  (2F10.5.4I2) 

TPRNT  ( 3) -TPRNT ( 1 ) 

TPRNT (5) -TPRNT ( 1 ) 

TPRNT(4)«2,*TPRNT(2> 

TPRNT(6)«2,*TPRNT(2> 

IF(J2)  45,45, 44 
READ(5,40)TARLP 
READ<5,40 )TA8LR 

READ  (5;40)  TO, PA. AR.XKN, ANR.TOPEN.EUSTAR.EPS 
READ (5, 40)  RbaR,RMIN.XNRAR,TFRAR|TSRAR,RHOS,GAM.CSTAR 
READ  <5;4Cn  ES. EG. CSCG.FEXT, ALPHA 
FORMAT  (8F10.5) 

TSOL -ALPHA/ (RRAR«RBAR) 

DO  30  1.1,6 

TPRNT(I)«TPRNT(!)/TSOL 

WRITE  (6,42)  TO.RRAR,ES,PA.RM!N,EG,AB,XNBAR,CSCG,XKN,TFRAR,FEXT. 

anr.tsbar. alpha, topen.rhos.elstar. GAM  .eps.cstar 

FORmaT(1H0,10X, ' INPUT  DATA*// 

e,Fl0,3,10x,6HRBAR  « , F10 . 3, 1®X , 6HES  «,F10,3/ 
»,F10f 3.1UX.6HRMIN  « , F10 , 3, 10X , 6HEG  «,F10,3/ 
-.F10.3.10X, 6HNBAR  »,F10.3,10X,6HCSCG  «,F10,3/ 
8,F1P,3,1HX,6HTFBAR«,F10,3,10X.6HFEXT  -,F1P,3/ 

_ _ _  *,F10I3I10X,6HTSBAR«,F10.3,10X.6HALPHA«.F1P,6/ 

20X.6HTOPEN-.F10 ,3, 10 X.6HRHOS  «,F10,3/ 
20X.6HLSTAR-,F10I3,10X,6HGAMMA-.F10.3/ 

20X.6MEPS  *,F10,3,10X,6HCSTAR»,F10.3) 

WRITE  (4,80) 

FORmaT-(1H0.3X, *  TIME* ,4X, ’PRESSURE* ,3X. *DLPDT*.7x, t RATE’ ,5X, 

*TFMPF*  ,5X.  »T£MPC*  ,5X»  •  TEMPS  • ,  6X ,  *  TL  AST  • ,  4X ,  *  DEPTH  • , 
».4X, ’EXPONENT*) 

CALL  EVFS  <N, TPRNT) 

I F ( J3 )  36,36,34 
IF(RATE>35,35.36 
ANR-0 ,5-ANR 
N-12 

GO  TO  41 

IF  (Jl)  60,60,50 
WRITE  <6,16)NAME 

CALL  XYPLOT  (J, 50, 100, 3, TIME, PRES, TEMP, BRATE) 

PRECEDING  PAGE  BUNK 


20X.6MTO 

20X.6HPA 

20X.6HAR 

20X.6HKN 

20X.6HANR 


L9,r«i0 

COATES 

FORTRAN  DECK  'MAIN 

i 

40 

60  GO  TO  IP 

41 

RETURN 

42 

END 

MESSAGE'  FOR  ABOVE  COMPILATION, 
VERSION  4  HOD  1 


p 


026 


10/24/69 


013,T«10  COATES 


FORTRAN  DECK  »SETUP 


i 


10/24/6 


SUBROUTINE  SETUP  <T,Y,SIG,N> 

DIMENSION  T<2),Y<20),SIG<20> 

COMMON  TPRNT<6),J1,J2,J3,NPTS.I, 

TO,PA,AB,XKN,ANR,TOPEN,ELSTAR,EPS, 

PRAR,Rm1N,XNRAR,TFRAR,TRRAR,RHOS,GAM,CSTAR, 

FS, EG, CSCG,FEXT, ALPHA, 

CAPGAm,PCRIT,GAMF1»GAMF2,GAMF3,TCHAM,TS0L,BB,DE,TCI» 
v, R,R,nLPDT, RATE, PRESS, TF,TC,TS,TLAST. DEPTH, AN, XNP, 
TARLP(4),TABLR(8>,RP 

COMMON  J,  TIME  1 1P>0) , PRES  1100) ,  TEMP  (1R0), BRATFI 100) 

REFERENCE  PRESSURE  IS  1000PSI 
REFERENCE  RATE  IS  RATE  AT  1000PSI 
REFERENCE  CONDITIONING  TEMP  IS  25C 
CALCULATE  CONSTANTS 

CAPGAMcGAM*(2,/(GAM+l«))**((CAM+l')/(GAM*l«) ) 
PCRTT«<?,/(GaM*1 . ))*•( GAM/ < GAM-1 , > ) 
GAMri«SORT{2,*GAM/<CAPGAM«(  GAM-1, ) > ) 

GAMF2«2./GAM 
GAMF3«(GAM+1. )/GAM 
AA»1.5*(1.0-XNRAR) 

BRsAA/(t  ,0-Aa ) 

DD«PHOS«XKN*CSTAR*RRAR/(1000,*32,17) 

TCHAM»ELSTAR/(CAPGAM«CSTAR*12.*DD) 

Cl*FG*1000,/(2,*l,9fl) 

in 

M*0 

T( 1 ) «0 , 

CALCULATE  APPROXIMATE  FLAME  TEMP 

T0-T0+273. 

TFRAD«FFXT*CSCG*TO 
TF«TFBAR*CSCG*<T0-298.)*TFRAU 
I F ( J2 )  *51,50,52 

CALCULATE  INITIAL  PRESSURE,  SUMmERFIELD  LAW 

xnp-xnrar 

RT»(TF/TFRAR)«*(1.*XNP)*EXP(-C1«C1,/TF-1,/TFBAR) ) 

P«( ( (1.0*RB)*RT*DO-1,)/RB)**1.5 
CC»PR*P*»0 , 667 
XNP«1,0-0,667*CC/(1.*CC) 

RP»P#<1,*RR)/(1.*CC) 

R«RT*RP 

ITFRATE  ON  TF 

tfrad»tfrad/r 

M«M*1 

IF (M-10)  26,26,27 
GO  TO  10 
GO  TO  60 

CALCULATE  INITIAL  PRESSURE,  VIELLE  LAW 

xnp.xnbar 

RT«<TF/TFPAR)**(1,+XNP)*EXP(-C1*(1,/TF-1,/TFBAR>) 

P«RT*DD** ( 1, /( 1 • -XNP) ) 

R«P»*XNP 
GO  TO  61 

CALCULATE  INITIAL  PRESSURE,  P,R  TABLE 

PclP00, 

r«rpar/pd 

CALI  PRATF(P,OP,XNP,TA0LP,TABLR,NPTS, I ) 
RT«(TF/TFRAR)**<1.+XNP)«EXP(-C1*(1./TF-1,/TFBAR) ) 

RP«RT*RP 


1 


19,Tal0  COATES  FORTRAN  DECK  'SETUP  » 


44 

IF( aRS(O«RP>«.0006)  35,36.30 

45 

30 

RaRP 

46 

P*1P0O.«RP*DD/RRAR 

47 

GO  TO  60 

46 

35 

P«P/1000, 

49 

RaR/RRAP 

50 

GO  TO  61 

C 

CALCULATE  SOL  10  TEMPERATURES 

51 

60 

TSBAR«TS.BAR*273. 

92 

TSaTSPAP/(1..1.98*TSBAR*ALOGCR>/(ESal000.)) 

53 

V<l>aTS/TO 

54 

Ttal ,13393«Y(1 13393 

55 

V(2)»?.*Y(1)*T1 

56 

DO  ?0  1*3.7 

57 

20 

Y(I)a?.*Y(I-l)»Y(I-?> 

56 

Y(8)«(VI7)*l.)/2. 

59 

Y(9»aP 

60 

Y(l0)aR 

61 

YC11 )«TF/TFBAR 

62 

Y(lP)aV(H) 

63 

TCI *TF 

64 

TCaTF 

65 

RETURN 

66 

END 

MESSAGES  FOR  ABOVE  COMPILATION. 


026 


10/24/6V 


C19,T«10  COATES 


FORTRAN  RECK  ‘PRATE 


i 


10/24/ftV 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


subroutine  prate <p,rp,xnp,tablp*tablr*npts* i > 

DIMENSION  TADLP(8»*  TAOl_R<6> 

1  IE  mBLP<I>*P>  2*4*4 

2  1*1-1 

IE  <  I  *1  >  3*  3# 1 

3  I  «1 

GO  TO  7 

4  IE  (TAHLPI » -P )  8*5*5 

5  !«U1 

IE  ( I-NRTS)10*6.6 

6  I *NP  TS*1 
GO  TO  7 

IP  IE  (TABt P<I*1>*P>  7,5*5 

7  XNP«ALOO(TARlR(  I*1)/TABLR<I>  >/AUOG(TABLP(  i*i>/tablpu  ) ) 

4  RP«TARLRC I )*(P/TABLP< I ) )**XNP 
RETURN 
END 


3  MESSAGES  EOR  ABOVE  COMPILATION, 


15»T«10  COATES 


FORTRAN  DECK  'MFEQ  « 


1 

2 

3 

3 

3 

3 

3 

3 

3 

4 

5 

6 
7 
6 
9 


10 

11 

C 

12 

13 

38 

14 

15 

40 

16 

C 

41 

17 

18 

10 

19 

20 

21 

22 

C 

20 

23 

C 

24 

25 

C 

26 

27 

C 

28 

29 

60 

30 

31 

48 

32 

33 

40 

34 

35 

50 

36 

37 

51 

38 

39 

55 

40 

C 

70 

41 

42 

C 

• 

43 

74 

44 

SMBPOUflNF  OfFEO  <T,Y.DYI)X,N,TPR) 

D ! MFNS I  ON  T<2)»Y(20),DYnX(20) 

COMMON  TPRNT(6),J1, J2,J3,NPTS,1, 

•  TO# PA, AB.XKN, ANR,TOPEN» ELSTAR.EPS, 

•  PRAR,RmIN,XNRaR,TFBAR,TSRAR,RHOS,GAM,CSTAR, 

•  rs. EG, rscc#FEXT, ALPHA, 

•  OApGAs,PCR ! T , GAMFl,GAMF2,GAMF3,TCHAH, TS0L,6B, DC, TCI  * 

•  X,P,R,OLpDT, RATE, PRESS, TF,TC,TS,TLAST,OEPTH, AN, XNP, 

•  TARLP(0),TABLR<8),RP 

common  J,TIME(lP»0),PRES<100),TEMP(100),BRATrC100) 

X«T(l)*TSOL 

P*  Y  ( 9 ) 

R«Y(10) 

RATF»R84R»Y<1#I) 

PRESS«100CI.#y<9) 

TF*TFRAR*Y<11) 

TC«TFRAQ*Y<12) 

CALCULATE  NOZZLE  AREA 
IF  # X«TOPFN)  30,30,40 
AN«AB/XYN*(ltR*(ANR-1.0)*X/TOpEN) 

GO  TO  41 

an*ab/xkn*anr 

A*AN/(On«A8/XKN) 

CHECK  FOR  EXTINGUISHMENT 
IF  (RATF*RMtN)  10,10,20 
RATr«0, 

N»0 

CALL  PRlNT(T,Y.nYOX,N,TPR) 

RETURN 
CONTI NUF 

CALCULATE  SURFACE  FLUX 

F$«P*mi)«l.*CSCGMTFBAR*T0-298,*TF>/T0>*FFXT 
CALCULATE  SURFACE  TEMP  CHANGE 
T1>Y(2)*0,267A6«FS/R 
DYDy<l)«56,»H*H*(Tl-2.0*Y(l)»Y(?) ) 

CALCULATE  RATE  CHANGE 
C2al000.*FS/(l,90«TO) 

DR»P*C2«DYDX(1)/<Y<1)*V<1)) 

CALCULATE  PRESSURE  CHANGE 
RNC«SORT(TF/TC!>*P»A 
PRATaPA/PRESS 

if<prat-pcrit)  55.55,4a 

PRATFaPRAT**GAMF2*PRATa*GAMF3 

IF  (PRATF)  49,49,80 

RSar, 

GO  TO  51 

H8«ONC*rPSaGAMFt*S0RT< PRATF) 

IF  fRS-PNC)  51,51,55 

RN«RS 

GO  TO  70 

RNaPNC 

DPaTSOL/TCMAM  *G AMaTC/TC I *H* < TF/TC*RN/R ) 

calculate  chamber  TEMP  CHANGE 

OTC» ( TSOL/TCham ) #TC/P*TC/TC I*R*(GAM*TF/TC*1.t(1,^GAM)«RN/R) 
IF<J2)  75,76,77 

CALCULATE  RP.XNP  SUMMERFIELD  LAW 
CCaop«P*»0,667 
XNP«l,0-Of667*CC/(l.*CC) 


10/24/69 


01*»T>10  corns 


FORTRAN  PECK  » D IFEO 


* 


10/24/69 


45 

RR»P*<1.*BR)/(1.*CC) 

46 

GO  TO  80 

C 

CALCULATE  RP * XNP  VIELLE  LAW 

47 

75 

XNPaXNBAR 

•  48 

RP»P**XNP 

49 

GO  TO  80 

C 

CALCULATE  RP# XNP  R# P  TABLE 

80 

77 

Palfl0O.*P 

?1 

CALI  PRATE (P,RP. XNP #TA8LP»TABLR»NPTS, 1 > 

52 

RPaPP/RPAP 

53 

P»P/1000« 

54 

GO  TO  60 

C 

CALCULATE  FLAME  TEMP  CHANGE 

55 

8P 

ETAb1.*VNP*10P0.*EG/<3,96*TF> 

56 

DTFaTF/FTA*(DR/R*XNP*DP/P) 

C 

calculate  solid  temp  change 

57 

DTDV(2)«42,*R*R*(Y(1 )-2.*V(2)*V(3) )*0.4615*(Y(3)-Y(1))*DR/R 

58 

DVDv<3)a30.*R*R»(V(?)-2.«Y(3)*Y(4))*0.d442*<Y(4)-.Y<2))*DR/K 

59 

DYDX(4)a20.«R«R*(Y(3)-2.*Y(4)*Y(5))-l.l323*{Y(5)-Y(3))*DR/H 

60 

DYDX(5)al2.*R*R«(Y(4)»2.*Y<5)4Y(6M*1.3020*(Y(6)^Y(4))*DR/R 

61 

DYDX ( 6 ) a  6,*R#R#(Y(5)-2.#Y(6)*Y(7M-1.3114*(Y(7)^Y(5n#np/R 

62 

DYDX<7)a  2,*R*R*(v(6)»2.*Y(7)*Y(8))*1.0619*<Y(8)*Y(6))*OR/R 

63 

DYDX(6)atP9l*R#R*(Y(7)-2**VC8>*l,)  •0.2357«<1,0-Y<7> )*DR/R 

64 

DYDX(9)«DP 

65 

DYDX(101iDR 

66 

DVOY<ll)«OTF/TFPAR 

67 

0YDXC12)«DTC/TFRAR 

68 

RETURN 

69 

END 

*0  MESSAGES  FOR  ABOVF  COMPILATION. 
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026 
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COATES 


FORTRAN  DECK  'PRINT 


t 
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1 

2 

3 

3 

3 

3 

3 

3 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 


SUBROUTINE  PRINT  <T,Y,DYDX#N,TPR> 

DIMENSION  T(2)#V(?0),DYnxi20),TPH<2) 

COMMON  TPRNT ( A ) , Jl, J2* J3,NPTS. I » 

•  TO, P A, aB.X*N,ANH,TOPEN, ELSTAB.EPS* 

•  obar,rmin,xnrap,tfbar,tsrar,rhos,gam,cstar, 

•  FS. EG, CSCG,FEXT, ALPHA, 

•  CAPGAh,PCRIT, GAMF1«GAMF2, GAMF3*  TCHAM, TSOL,bB,DC*TC I , 

•  X,  p,  R, nLPDT, RATE, PRESS »TF,TC#TS,TL AST .DEPTH, AN* XNP, 

•  TARIPC A>*TA8LR<8),RP 

COMMON  j. T 1  Me U00 >, PRES < 100), TEMP* 100 >,BR ATE (100) 

TS«T0*Va>-273. 

TLAST»TO*Y<8)-273. 

DEPTH* ALPHA/RRAR*2.593/R 
DLPnT»DVDX(9)/(P*TS0L) 

SRATE*RP*RBAR 

14  WRITE  (A, 15)  X. PRESS, DLPDT#RaTE.SRaTE,TF,TC.TS,1iaST,CEPTH, AN, XNP 

15  FORMAT  CF10.5.F10.2,3F10,4,2F10.1,2F10.2,2F10,S,F10.4> 
TSTop«TprnT(2)*TS0L 

IF  (X.GF.TSTOP)  N»0 

J*J*1 

TIMF< J)«x 

PRESC J)«P 

TEMPI J)«Y(12) 

0RATEI JI«R 

RETURN 

END 
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The  objective  of  this  project  was  to  carry  out  both  experimental  and 
analytical  studies  leading  to  improved  design  procedures  for  predicting  thrust 
termination  of  single-chamber  controllable  solid  motors.  The  experimental  work 
consisted  of  measuring  quantitatively  the  effects  of  incident  thermal  radiation  on 
low  pressure  burning  rates  and  deflagration  limits  of  typical  solid  propellants. 

A  technique  was  developed  in  which  small  cylindrical  samples  of  the  propellant 
were  burned  inside  an  electrically  heated  tube  furnace.  At  furnace  wall  tempera-^ 
tures  corresponding  to  the  interior  surfaces  of  a  rocket  motor,  the  burning  rates 
were  observed  to  be  as  much  as  50%  greater  at  lOpsia  than  they  were  when  measured 
in  the  conventional  strand  burner  configuration.  Analytical  work  was  performed 
using  an  improved  mathematical  model  of  the  transient  combustion  process  of  a 
solid  propellant. (  This  study  led  to  the  conclusion  that  experimentally  character¬ 
izing  the  extinguifcptoility  of  a  propellant  in  terms  of  a  critical  dp/dt  provided 
very  little  useful  design  information.  On  the  other  hand,  characterizing  exting¬ 
uishing  in  terms  of  thii  product  L*r  2  was  shown  to  provide  a  design  criteria  with 
the  promise  of  very  general  applicability. 


X 


Solid  Propellant 
Cori>ustion 
Rocket  Motor 
Characteristic  Length 
Burning  Rate 
Extinquishaent 
Depressurization 
Depressurization 
Deflagration  Liait 
Nozzle  Area 


